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OCEAN  BOTTOM  BREAKOUT  FORCES 


Including  Field  Test  Data  and  the  Development  of  an  Analytical  Method 
Technical  Report  R- 591 
by 

Bruce  J,  Muga 


ABSTRACT 

Theoretical  and  experimental  studies  were  conducted  in  order  to 
arrive  at  some  appropriate  engineering  estimates  of  the  force  required  to 
extract  bodies  of  various  sizes  and  shapes  from  the  ocean  bottom  sediments, 
A  review  of  the  literature  concerned  with  breakout  forces  Is  presented, 

From  tests  conducted  In  San  Francisco  Bay  on  variously  shaped 
objects  having  submerged  weights  of  up  to  22,200  pounds,  it  was  found  that 
the  following  empirical  formula  described  well  the  breakout  force  require¬ 
ment: 


F 


0,20  AmiXqde 


-0,00840  (t  -  2801 


where 


F  *»  breakout  force,  lb 

Ami>  “  hurlzortal  projection  of  the  maximum  contact  urea,  In,* 

qd  »  average  supporting  pressure  provided  by  the  soil  to  maintain 
the  embedded  object  In  static  equilibrium,  lb/ In.2 

t  ■  time  allowed  for  breakout,  min 


'.umciw  .  .  .  Results  yielded  by  the  empirical  formula  agreed  very  well  with  a 

jnwn  *,nn  '■* W  molested  theoretical  procedure  bused  on  an  Iterative  solution  of  a  lumped 
"  •fwfmriater  modol  of  the  ocean  bottom,  Solutions  worn  obtained  for  various 
I ld#d  conditions  and  bottom  (ob|oct)  geometries. 
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INTRODUCTION 

There  Is  an  increasing  interest  being  shown  In  deep  sea  operations. 

This  Interest,  having  scientific  origins,  has  been  uugmarited  by  the  need  to 
enhance  the  rational  security  and  to  prevent  catastrophes.  As  a  result,  the 
engineering  aspects  of  certain  practical  problems  of  Immediate  Interest  have 
tended  to  dwarf  the  longer  terrr  scientific  aspects,  The  study  reported 
herein,  which  was  set  Into  motion  following  the  loss  of  the  USS  THRESHER, 
is  concerned  with  one  of  these  practical  problems. 

Briefly,  the  objectives  of  the  study,  appropriately  described  as  the 
"breakout  force"  problem,  are  threefold.  First,  some  appropriate  engineering 
estimates  must  be  made  of  the  force  required  to  extract  bodies  of  various 
shapes,  sizes,  and  skin  surface  compositions  from  the  ocean  bottom  sedi¬ 
ments,  Second,  it  must  be  determined  how  far  the  estimates  may  be  In  error. 
Third,  the  effectiveness  of  various  schemes  for  reducing  the  required  force 
must  be  evaluated. 

The  study  Is  a  continuing  3-year  effort  of  which  this  report  ia  intended 
to  describe  the  progress  made  In  the  first  2  years.  Thus,  although  this  report 
la  not  considered  to  be  a  complete  treatment  of  the  subject  matter,  It  makes 
It  possible  to  estimate  the  force  requirements  with  some  degree  of  confidence. 
In  addition,  the  theoretical  treatment  has  been  substantially  verified  by  field 
tests  and  Is  believed  to  be  a  significant  contribution  to  analytical  methods 
which  heretofore  have  not  been  widely  used,  It  appears  that  many  practical 
problems  exist  for  which  such  methods  have  application, 


BACKGROUND 

Documented  experience  relative  to  the  problem  of  freeing  objects 
from  the  ocean  bottom  Is  limited.  The  reports  referenced  In  Table  1  mention 
the  "mud  suction"  force;  however,  no  measurements  of  such  a  force  seem  to 
have  been  made.  A  description  of  the  nature  of  the  problem  may  be  pro¬ 
vided  by  brief  reviews  of  the  salvage  operations  listed  in  Table  1. 
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Table  1.  Selected  Marine  Salvage  Operations 


Sunken  Objects 

Date  of 

Salvage  Operations 

Reference 

SODRA  SVERIGE 

1896-97 

Anonymous  (1956) 

LIBERTE 

1911-25 

J.  Crapaud  (1925) 

S-51 

1925-26 

U.  S.  Navy  (1927) 

S-4 

1927-28 

U.  S,  Navy  (1929) 

USS  SOUALUS 

1939 

(a)  R,  A.  Tusler  (1940) 

(b)  C,  E,  Momsen  (1964) 

USS  LAFAYETTE 
(ex  SS  NORMANDIE) 

1942-43 

U.S,  Navy  (1946) 

PHOENIXES 

1944 

E.  Ellsberg  (1960) 

SODRA  SVERIGE.  ThoSODHA  SVERIGE,  a  cargo-passenger 
steamship  with  an  ROC-ton  displacement,  stink  In  the  Baltic:  Sea  in  1895  In  a 
depth  c*  185  feet.  The  ship  came  to  rest  at  a  sharp  angle  from  the  vertical 
and  during  the  course  of  a  year  sank  about  10  feet  Into  the  day  bottom, 
Calculations  Indicated  that  the  ship  had  a  submerged  weight  of  600  tons 
and  that  a  force  of  960  tons  would  be  sufficient  to  break  It  loose  from  the 
bottom.  Sixteen  wooden  pontoons,  each  having  a  lifting  force  of  60  tons, 
were  attached  to  the  ship  and  pumped  out.  This  was  sufficient  to  rlghten 
the  ship  and  raise  It  off  the  bottom, 

LIBERIE.  The  French  battleship  LIBERTE,  with  an  B, 000-ton 
displacement,  sank  In  the  harbor  at  Toulon  In  191 1 ,  The  selvage  effort 
extended  over  a  period  of  14  years,  During  the  long  period  of  submergence, 
the  wreck  settled  Into  the  mud  and  Crapaud  (1925)  reports  that  "a  con¬ 
siderable  part  of  the  task  of  the  salvors  consisted  In  breaking  this  contact 
and  freeing  the  bulk  so  that  It  could  be  lifted  and  towed  away,"  No 
Information  was  given  to  permit  an  estimation  of  the  breakout  force. 

S-61.  The  S-51  wasa  1,230-ton  submarine  which  sank  In  1926 
approximately  14  miles  east  of  Block  Island,  The  depth  at  the  cite  was 
132  feet  and  the  submerged  weight  was  estimated  to  be  1,000  tone,  The 
boat  came  to  rest  on  a  clay  bottom  with  on  1 1-degree  port  list,  Ellsberg 
( 1927)  estimated  that  the  breakout  force  was  about  8,000  tone,  "ii  force  so 
large  we  could  never  hope  to  overcome  It  by  direct  lift,"  His  plan,  which 
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vA/ac  ovominH  <;i irrAcqfi ii ly  wfis  tn  "break  the  suction  bv  lettinq  water  in 
between  hull  and  clay  in  two  ways— first,  by  rolling  the  boat  to  starboard, 
and  second  by  lifting  her  one  end  (stern)  first," 

S-4.  The  S-4  was  an  830-ton  boat  which  sank  in  1927  in  102  feet  of 
water  off  Provincetown.  It  was  initially  and  intermittently  buried  in  a  very 
permeable  mud  to  a  depth  of  7  or  8  feet.  The  boat,  whose  submerged  weight 
amounted  to  722  tons,  was  raised  by  lifting  the  stern  first,  Saunders  ( 1929) 
states  that  "the  bottom  had  an  upper  layer  of  very  soft  silt  or  mud  not  more 
than  one  foot  deep.  Underneath  this,  ths  bottom  was  more  soft  than  hard, 
of  a  decidedly  sand  character,  mixed  with  minute  shells.  The  texture  of  the 
bottom  was  sufficiently  coarse  to  permit  the  passage  of  water  through  It. ..yet 
sufficiently  firm  to  hold  it3  position  when  excavating  tunnels  underneath  the 
vessel.  Due  to  the  permeable  characteristic  of  the  bottom  It  Is  estimated  that 
the  so-called  'suction  effect'  on  the  S-4  was  practically  nil,"  As  a  matter  of 
fact,  there  are  no  indications  that  breakout  was  a  problem, 

USS  SQUALUS.  The  salvage  of  the  USS  SQUALUS  Is  perhaps  the 
most  widely  reported  and  documented  of  all  marine  salvage  operations,  The 
USS  SQUAI.US,  which  sank  In  1939  about  5  miles  south  of  the  Isles  of  Shoals 
off  Portsmouth  Harbor,  was  a  1 ,450-ton-dlsplacement  boat  having  a  sub¬ 
merged  weight  of  1,100  tons.  The  boat  came  to  rest  with  a  10-dogree-up 
angle  In  240  feet  of  water  on  a  mud  botiom  In  which  the  stern  was  burled 
up  to  the  superstructure  deck,  The  entire  operation  consisted  of  five  separata 
attempted  lifts,  three  of  which  were  from  a  mud  bottom,  Only  tho  first  lift 
Is  pertinent  to  this  report,  No  estimate  of  the  breakout  force  Is  reported; 
however,  in  a  review  of  the  events  describing  the  unsuccessful  lift  of  13  July, 
Tuslcr  (1940)  Indicates  that  the  “unknown  amount  of  mud  suction  tending 
to  hold  the  bow  down"  was  one  of  the  main  factors  contributing  to  the 
failure,  Previous  to  the  attempted  lift,  Tusler  states  that  "the  bow  had  sunk 
down  an  unknown  amount  Into  the  soft  mud  of  the  bottom,  hut  due  to  the 
shape  of  the  bow,  It  was  thought  that  the  mud  suotlon  would  be  relatively 
Insignificant,"  In  any  event  during  each  attempted  lift,  the  USS  SQUALUS 
was  raised  by  lifting  one  and  first, 

USS  LAFAYETTE  (ex  SS  NORMANDIE),  The  SS  NORMANDIE 
was  a  06,000-ton  passenger  vossul  which  sank  In  49  foot  of  water  adjacent 
to  Pier  88,  Mew  York  City  Harbor,  The  submerged  weight  was  estimated  to 
bo  60,000  tons,  The  vossol  come  to  rust  lying  on  one  uldu  In  an  organic  river 
mud  which  woa  ubout  26  foot  thick  and  which  was  underlain  by  «  gray 
organic  silty  cloy  having  a  compressive  strength  of  from  0,3  to  0.6  ton/ft* , 
(This  operation  Is  of  Interest  since  It  appears  to  have  boon  tho  first  tlmo  tho 
principles  of  soil  mechanics  wore  considered  In  u  salvage  oporatlon  of  this 


type,  i'he  fact  that  Professor  K.  Terzaghi  served  as  consultant  is  a  historical 
note  of  some  Importance.  Ho  was  not  concerned  with  breakout  but  rather 
with  determining  If  the  soil  had  sufficient  strength  to  support  the  vessel 
without  serious  movement  during  the  year  and  a  half  of  preparation  for 
righting,  with  the  effect  dredging  would  have  on  the  settlement  of  the  ship, 
and  with  the  probable  effective  soil  bearing  plane. )  It  was  anticipated  that 
breaking  contact  between  the  ship's  hull  and  the  mud  would  be  a  serious 
problem.  Accordingly.  In  addition  to  pumping  out  some  16,000  tons  of 
mud  which  had  entered  the  hull  through  the  cargo  doors  and  portholes, 
numerous  porthole  patches  were  fitted  with  pipes  through  which  water  and 
compressed  air  could  be  jetted  to  disintegrate  tho  mud.  The  flotation  of  the 
vessel  was  preceded  by  a  rotation  or  turning  operation,  Masters  (1964)  notes 
that  during  the  rotation  operation,  the  air  and  water  lots  were  sot  to  work, 
although  the  vessel  did  not  stick  bs  expected. 

PHOENIXES.  The  PHOENIXES  were  200-foot-long  floating  blocks 
of  reinforced  concrete  which  were  to  be  sunk  In  a  line  off  the  Normandy 
beaches  to  provide  a  breakwater  during  the  European  Invasion,  Each  unit 
was  60  feet  wide,  60  feet  high,  and  displaced  6,000  tons,  They  were  divided 
Into  watertight  compartments  and  fitted  with  valves  tor  controllable  flooding. 
Approximately  100  of  these  PHOENIXES  were  purposely  sunk  In  staging 
areas  off  the  south  coast  of  England  prior  to  Invasion,  Tha  first  attempt  to 
ref  lost  e  PHOENIX  by  pumping  felled,  It  was  determined  that  tha  mud 
bottom  euotlon  wee  holding  tha  PHOENIX  down,  The  traditional  method  of 
breaking  thle  contaot  Is  to  apply  buoyancy  m  one  ond  and  to  use  the  ship  as 
•  lever,  In  thle  manner  the  contact  la  brokan  along  tha  bottom,  eventually 
freeing  tha  vessel  from  the  mud,  However,  In  tha  case  of  the  PHOENIXES 
there  was  not  enough  time  to  allow  this  system  to  work,  because  the  rising 
tldas  would  submargu  tha  pump  platforms  before  the  lever  action  could  be 
made  affective,  A  seoond  alternative  Is  to  |et  air  or  wator  underneath  the 
sunken  vessel  to  partially  reduco  the  contact  end  lessen  tha  holding  foroa,  In 
tha  oase  of  tha  PHOENIXES,  compressed  air  was  employed  (all  pumps  wore 
being  used  to  empty  the  flooded  compartments)  to  reduce  the  holding  force 
to  allow  the  excess  buoyancy  to  float  the  PHOENIXES, 

In  addition  to  tha  case  histories  cited  above,  thorn  are  many  other 
records  of  maritime  salvage  operations  In  which  a  ship  has  beun  raised  from 
a  mud  bottom  under  very  unfavorable  circumstances,  The  background 
provided  heroin  Is  nut  Intended  tv)  bo  an  all-inclusive  treatment  of  salvage, 
Only  those  oases  In  which  the  breakout  foroa  was  alluded  to  In  the  published 
literature  wore  selected,  It  Is  worthy  of  note,  however,  that  tho  published 
record  (Bowman,  19t34)  of  the  salvage  of  tho  entire  Qermon  High  Seas  fleet, 
which  woe  scuttled  at  Scapa  Flow,  does  riot  mention  that  broukau1  was  a 


problem.  Further,  the  U.  S.  Navy's  experience  with  the  ex-German 
submarine  1 105  did  not  disclose  the  breakout  problem,  although  the  tests 
were  carried  out  with  the  submarine  eventually  lying  on  a  mud  bottom 
(lotter  from  Commander  Task  Unit  49,4,3), 

In  summary,  the  problem  of  freeing  vessels  !.'w  *ho  ocern  bottom 
has  been  recognized  for  a  long  time,  However,  a  number  of  factoi  s  In 
combination  with  each  other  have  tended  to  prevent'the  synthesis  of  the 
basic  principles  of  naval  architecture  and  those  of  soil  mechanics.  The 
development  of  modern  soil  mechanics  dates  from  the  work  of  K,,  Terz^ghl 
during  the  1920s  and  Is  a  relatively  recent  event  In  torim  ol  men's  experience 
with  sea-going  vessels,  Moreover,  when  the  unknown  behavior  of  ocean 
bottom  sediments  Is  considered,  one  can  sppreclste  why  the  marine  salvor 
attaches  more  Importance  to  hydrostatic  calculations  than  to  gross  estimates 
of  the  mud  suotlon  affect, 


FIELD  TESTS 

A  field  vest  program  wai  deilgnod  to  correlate  the  breekout  force  with 
the  breakout  time,  object  else,  object  shape,  end  soil  strength.  The  breakout 
force  calculated  from  the  resulting  empirical  formula  was  compared  with  the 
analytical  value.  A  detailed  description  of  the  site  selection  end  soil  borings 
Is  given  In  Appendix  A,  A  detailed  description  of  the  selection  of  boOy  shepee 
end  dimensions  Is  given  In  Appendix  B, 

Limitations 

A  number  of  faatori  disclosed  by  «  preliminary  survey  end  recant 
experience  Imposed  some  limitations  on  the  conduct  of  the  field  tests,  These 
factors  were: 

1,  The  tests  could  not  be  conducted  near  -any  existing  pier,  wharf,  or  other 
fixed  structure  due  to  the  presence  of  eheh  fragments  end  other  debris  In  the 
soil, 

2,  Each  tost  was  required  to  ba  conducted  Ir.  an  essentially  undisturbed 
deposit,  even  though  local  remolding  of  the  soli  was  known  to  take  place, 

3,  Since  breakout  Is  essentially  *  time-dependent  process,  It  Is  necessary  to 
apply  relatively  constant  force  levels  over  long  periods  end  to  measure  the 
movement  of  the  ob|«ot  with  respect  to  the  undisturbed  deposit  In  v.l  luh  it 
Is  embedded, 


f  i 

; 

i  ( 


Apparatus 

The  use  of  a  mobile  bottom-resting  platform,  wen  if  one  had  been 
available,  was  determined  to  be  uneconomical.  A  proposal  to  drive  groups  of 
piles  at  various  locations  within  San  Francisco  Bay  also  proved  to  be 
uneconomical.  It  was  Intended  that  the  piles  would  support  fixed  decks 
from  which  the  tests  could  be  conducted. 


matar  tor  maaauUng  dlaplactmant 
(A^l  of  objtct  rilatlva  to  bare* 


"04fp«nt»r"  itop  . 
lotddynimomttir  ■ 


oountwwtlflht  load 


ct«<r 


moor'n«"_n* 


watar  aupply  ayitam  tor 
oounttfwalaht  loading  tanka 


'  motar  for  maaaurlng 
I  dlaplactmant  (Aj) 

|  of  barga  vtlatlva  to 
I  bottom 


A  •  i  A2,  whara  A  la  tha  dlaplacamant 
ol  tha  objaot  rolatlva  to  tha  oottnm 


Figure  1.  Schimitlo  of  ta»t  apparatua. 
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A  db-fout  x  MB- foot  floating  barge,  called  thy  NCEL  warping  tug,  was 
selected,  largely  because  of  its  nvniLbility.  The;  warping  '  1  ’  was  modified  to 
accommodate  a  counterweight  loading  system,  a  schematic  of  which  is  shown 
in  Figure  1  The  counterweight  system  permits  the  application  of  relatively 
constant  loads  over  long  periods  of  time,  since  it  compensates  for  tide  changes, 
varying  freeboard,  and  the  short-term  heave,  pitch,  and  roll  of  the  barge.  The 
Systran  consists  of  a  frame  composed  of  pontoons  (b  x  b  x  7  feet)  which  pro¬ 
vide  an  elevated  support  for  the  counterweight.  The  counterweight  proper 
consists  of  twin  (Moot-diameter  by  12-foot-long  tanks,  These  tanks  when 
filled  with  varying  amounts  of  water  tire  capable  of  providing  from  9,000  to 
<19,00(1  pounds  of  force.  The  force  levels  can  he  carefully  controlled,  A 
photograph  of  the  system  is  shown  in  Figure  2. 

The  "curpontnr  stopper, '  shown  in  Figure  3,  allows  an  adjustment  of 
tho  length  of  the  lino  from  the  counterweight  to  the  test  object,  It  is  a 
critical  piece  of  equipment  since  It  can  bo  used  to  corinoct  one  lino  anywhere 
along  tho  length  of  another  line  without  damaging  either  lino,  It  develops  the 
full  strength  of  tho  linos  and  does  not  disengage  when  the  lines  are  slockod, 


/ 


Figure  2,  Photograph  of  test  apparatus 
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i  he  warprng  lug  was  secured 
in  position  with  a  two-point  moor, 
fore;  and  aft.  The  stern  lino  was 
anchored  with  a  9, 000-pound  STATO* 
anchor  and  the  bow  line  was  made  fast 
to  tne  end  of  a  pier,  Sufficient  tension 
was  developed  in  the  mooring  lines  to 
resist  significant  transvorso  movements 
of  the  barge,  Position  was  checked  by 
measuring  (with  a  transit)  the  angle  of 
the  n  nor  Tintivn  to  the  pier  and  by 
metering  the  bow  line. 

Force  was  measured  and 
racordad  by  on  In-line  calibrated  strain- 
f  age  dynamometer.  Displacement  of 
the  test  specimen  relative  to  the  barge 
was  measured  by  a  potentiometer- 
backed  wheel  driven  ay  the  sheave 
over  which  the  lifting  line  was  reevod. 
Displacements  of  tho  barge  relative  to 
tho  bottom  wore  also  measured  by  o 
pov.ontlometor-backed  wheel  over 
which  a  small  counterweight  and  a 
mud  float  were  suspended  (see  Figure 
1),  Some  difficulty  was  experienced 
from  floating  debris,  which  tended  to 
become  entangled  In  the  wire  line,  A 
facsimile  of  an  oscillogram  for  one  of 
the  tests  Is  shown  In  Figure  4. 
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RESULTS  OF  FIELD  TESTS 


Presentation  of  Data 


The  data  as  reduced  from  the 
oscillograms  are  portrayed  graphically 
in  Figures  5  through  10.  It  Is  to  be 


emphasized  that  these  are  essentially 
raw  data  which  have  not  been 
corrected  for  oovlous  errors,  although 
some  test  data  were  discarded  primarily 
for  reasons  of  Incompleteness,  The 
data  are  presented  In  the  form  of  depth 
of  embedment  In  Inches  versus  time  In 
minutes  from  the  beginning  of  oach 
Ust.  The  beginning  of  each  test  was 
taken  to  be  the  time  at  which  the 
counterweight  tanks  began  to  be  filled 
with  water,  For  each  of  the  tests, 
unless  so  Indicated,  the  rate  of  load 
application  (that  is,  the  filling  of  the 
tanks)  was  uniform,  although  It  varied 
from  test  to  test,  For  each  test,  the 
opplled  load  over  and  above  the 
submerged  welpht  of  the  test  specimen 
(In  aalt  water)  is  Indicated,  In  addition, 
certain  other  Information  Is  presented 
for  each  test,  This  Includes  (1)  the 
time  elapsed  from  the  placement  of  the  objeot  to  the  start  of  the  test;  (2)  the 
location  of  the  test  object  relative  to  the  pier  In  bearing  and  distance;  and 
(3)  thu  local  shear  strength  at  o  depth  equivalent  to  the  depth  of  settlement 
of  the  object.  The  latter  was  obtained  by  divers  using  a  vane  device.  To  aid 
In  Interpreting  thu  figures,  the  following  example  Is  described  In  detail. 

Example  (Teet  No,  38,  Figure  6),  The  test  object  was  a  cube  which 
was  placed  In  position  240  fact  from  the  pier  at  an  angle  of  60  degrees 
measured  olookwlse  from  an  Imeglnar,  line  perpendicular  to  the  longitudinal 
axis  of  the  pier  and  having  Its  origin  at  the  southeast  corner  of  the  pier,  The 
object  settled  a  distance  of  36  Inches  Into  the  bottom  sediments  and  wbb 
allowed  to  remain  In  place  undisturbed  for  68  hours,  In-sltu  vane  shear 
strength  at  a  depth  of  36  Inches  was  measured  to  bn  0,83  psl.  The  submerged 
weight  (In  salt  water)  of  the  cube  was  20,000  pounds,  At  the  beginning  of 
the  test  a  load  of  9,000  pounds  was  applied  to  the  object,  which  corresponds 
to  the  dead  weight  of  the  counterweight  tanks,  Water  was  added  at  about  the 
rate  of  216  Ib/mln,  After  sipproxlmately  100  minutes,  the  total  applied  force 
was  30,600  and  the  object  moved  noticeably  and  broke  free  of  tho  bottom, 
The  data  from  all  of  the  tests  are  conveniently  displayed  In  Table  2, 

To  bring  some  order  out  of  the  large  number  of  parameters  required  to 
describe  the  breakout  phenomenon,  the  original  raw  data  presented  In 
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Figure  4.  Fioilmlle  of  teit  date. 


Figures  5  through  8  have  been  replotted  and  appear  in  Figures  *  1  Li  u  uiftjn  )  I 
The  ordinate  has  been  normalized  by  dividing  the  instantaneous  embedded 
depth  by  the  original  embedded  depth,  It  was  expected  that  the  freeing  of 
the  object  or  alternatively  the  failure  of  the  soil  medium  could  be  clearly 
Identified  and  correlated  both  for  similar  and  dissimilar  objocts.  However, 
this  hope  did  not  materialize, 

Formulation  of  Empirical  Equation 

Unfortunutely,  the  various  classical  theories  of  failure  which  have 
been  widely  employed  in  the  practice  of  soil  mechanics  and  foundation 
engineering,  such  as  the  Mohr-Coulomb  theory,  assume  that  the  stress  condi¬ 
tions  alone  determine  the  state  of  failure  of  a  material,  irrespective  of  the 
load  duration  arid  the  stress  history,  However,  It  ha»  become  Increasingly 
clear  that  the  load  duration  is  a  major  factor  In  the  breakout  process.  Schmid 
and  Kltago  (1666)  stated: 

"That  clay  soils  have  time-dependent  sheer  properties  has  been 
recognized  for  a  long  time,  but  the  profession  still  frequently  approaches 
problems  of  time-dependent  stress-strain  behavior  of  clay  soils  as  If  they  could 
be  analysed  completely  on  the  basis  of  the  theory  of  elasticity  or  plasticity 
without  regard  to  the  actual  rheological  properties,  Time  affects  are  often 
discussed,  If  at  all,  merely  on  a  qualitative  basis,  One  reason  for  this  may  bo 
the  complexity  resulting  when  time  Is  considered  an  additional  variable  In 
any  problem.  Another  reason  may  be  the  difficulty  of  abandoning  old  and 
familiar  concopts  that  usually  are  sufficient  for  most  structure!  materials,  and 
another  the  lack  of  familiarity  of  most  soil  engineers  with  rheological  theory, 
since  Its  application  to  problems  of  soli  mechanics  Is  still  developing. 

"There  appears  to  be  general  agreement,  however,  that  saturated  clay 
soils  do  behave  like  viscoelastic  or  viscoplastic  materials,  As  a  consequence, 
the  classical  failure  theories, ..cannot  and  do  not  completely  describe  the 
material  behavior  of  clay  soils,  Either  they  have  to  be  modified  to  permit 
.i  quantitative  assessment  of  stress  history,  temporaturo,  and  rate  of  loading 
or  they  have  to  be  replaced  by  theories  that  Include  these  effects." 

In  the  classical  theories,  one  or  two  of  such  material  parameters  as 
yield  stress,  Young’s  modulus,  or  Poisson's  ratio  are  sufficient  to  describe 
the  behavior  of  an  isotropic  materiel.  Such  simplifications  are  Inadequate  to 
predict  breakout  behavior, 

To  keep  the  problem  as  simple  as  possible  and  yet  retain  ell  of  the 
essential  features,  the  breakout  test  data  will  be  analyzed  stepwise  In  the 
following  paragraphs, 
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Beginning  with  a  very  simple  formulation  of  Ihe  mechanics  of 
breakout,  one  may  set* 


F  -  k  C  A 


(1) 


where 


F  =  breakout  force 

C  -■=  cohesion,  or  alternatively  a  measure  of  the  vane  shear  strength 

A  =  horizontal  projection  of  the  contact  area 

k  -  constant  which  is  a  function  of  object  size,  object  shape,  time 
duration  of  applied  force,  rate  force  is  applied,  soil  sensitivity,** 
and  the  elapsod  time  which  the  object  has  been  in  place  after 

the  Initial  disturbance 


Thus,  we  may  write 


or 


(2) 


Letting  C  and  k  take  on  slightly  different  meanings,  we  may  write 


£A 


(3) 


whero  C  ■  effective  average  cohesion  along  the  failure  surface  at  the  instant 
of  breakout 


Q  ■  constant 

R  °  slope  of  the  "failure  line"  when  log  (F/£A)  Is  plotted  versus 
time,  t 

t  =  time  allowed  for  breakout,  or  alternatively  the  elaosed  time 
during  which  the  breakout  force  Is  applied 


t0  ■  reference  time  in  minutes 


*  The  reader  Is  referred  to  the  nomenclature  on  page  31 . 

**  Ratio  of  cohesion  of  undisturbed  soil  to  cohesion  of  disturbed  soil  at  constant  water 
content, 
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Tha  constants  Q  and  R  are  functions  of  the  load  duration  or  strain 
rate.  In  Equation  3,  when  t  -  »*,  the  force  required  for  breakout  is  a  mini¬ 
mum.  Conversely,  ast  is  allowed  to  approach  zero,  that  is,  as  the  time 
allowed  for  breakout  becomes  increasingly  short,  the  force  requirement 
reaches  a  maximum  constant  value. 

The  quantity  C  requires  some  comment,  since  it  is  also  a  time- 
dependent  function  which  is  related  to  the  soil  sensitivity.  It  may  be 
estimated  by  an  equation  of  the  type 


where  s  •  degree  of  soil  sensitivity 

b  ■  numerical  constant  used  to  force  C  *  C  for  very  large  t,  In 
keeping  with  our  knowledge  of  thixotropic  material  behavior 

t,  *•  reference  time  related  to  the  tnixotroplc  behavior  of  a  material 
in  regaining  a  stated  percentage  of  its  strength  after  initial  dis¬ 
turbance 

To  Illustrate  that.  Equation  4  is  approximately  correct,  we  note  that  (1)  for 
t  “  0,  then  £■  C/»,  (2)  for  t  ■  tp  Q"  C,  and  (3)  for  very  large  t,  C  ■  C  even 
for  relatively  small  values  of  b.  Dimensionless  graphs  of  Equation  4  are 
shown  In  Figures  15  and  16  for  values  of  b  ••  1,0  and  s  ■  5.0  and  8.0, 
respectively.  Experimental  Information  on  the  validity  of  Equation  4  seems 
to  be  nonexistent,  In  addition,  the  reference  timet!  seems  to  be  highly 
variable,  being  very  short  (that  Is,  on  the  order  of  minutes)  for  such 
thixotropic  materials  as  drilling  muds  and  perhaps  very  long  (that  Is,  measured 
by  geologic  time)  for  many  deep  marine  sediments, 

On  the  basis  of  experimental  test  results  wo  may  estimate  »,  C,  t, ,  t0, 
and  the  constants Q  and  R,  and  then  compute  the  force  F  required  to  extract 
the  specimen  os  a  function  of  time,  t.  For  example 

F  -  QAQen,,-V  (6) 

or  for  maximum  C 

F  -  CAQe-R<,-‘o>  (6) 

It  Is  to  be  emphasized  that  the  reference  times  are  those  determined  from 
large-scale  field  tests, 
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Figure  16.  Dlmamlonleu  graph  of  Equation  4  Figure  16,  DlmanBlonleu  graph  of  Equation  4 
for  b*  1,0,  a  ■  6,0,  for  b  “  1,0,  i  ■  8,0. 


The  foregoing  Illustrates  how  the  field  data  may  be  analyzed  and  used 
for  predicting  breakout,  assuming  that  scale  effect  Is  negligible,  In  connection 
with  the  data  reduction  it  was  found  that  the  coneslon,  C,  as  obtained  by  vane 
shear  tests,  showed  marked  variability.  Thus,  as  an  alternative  measure  of  the 
sediment  strength,  It  was  expedient  to  use  the  quantity  qd,  which  Is  defined  as 
the  average  supporting  pressure  provided  by  the  soil  to  maintain  the  embedded 
object  In  static  equilibrium,  In  all  of  the  data  reduction  since  It  exhibited  very 
consistent  trends,  In  a  sense  this  Is  fortunate  since  the  problem  then  becomes 
completely  determinate,  being  no  longer  dependent  on  external  measurements, 

Soma  of  the  data  summarized  In  Table  2  are  presented  In  Figures  17 
and  18,  Figures  18a,  18b,  and  18c  ore  semllogarlthmlc  graphs  of  the  data 
appearing  In  Figures  17a,  17b,  and  17c,  respectively,  The  elapsed  time  over 
which  the  maximum  force  was  applied  appears  as  the  ordinate  In  all  figures, 

Iri  Figures  17a  and  17b  the  abscissa  Is  the  dimensionless  quantity,  F/(AmiXqd); 
In  Figure  17c  the  abscissa  Is  F/(Aqd),  In  Figures  17b,  17c,  18b,  and  18c  the 
force  F  represents  the  net  breakout  force,  which  Is  the  applied  force  minus 
the  submerged  weight  of  the  object,  not  only  that  portion  submerged  in  salt 
water,  but  also  that  portion  of  the  object  embedded  In  the  bottom  sediment. 

Figure  1  Qa  Indicates  S(.,me  of  the  tronds  in  selected  data  from  Table  2. 
The  coefficients Q,  R,  and  t0  used  In  Equation  5  are  ulso  shown  on  the  figure. 
Although  the  data  aru  extremely  limited,  the  figure  indicates  that  the  forces 
required  to  extract  the  cube  and  the  prism  are  higher  than  those  for  the 
cylinder  and  sphere,  This  seems  to  be  In  agreement  with  previous  field 
experience, 


Figure  18b  exhibits  the  clearest  consistent  relationships  for  all  of  the 
included  data.  That  is,  the  trends  sliuwn  in  Figure  10a  for  the  various  shaped 
objects  appear  to  be  almost  entirely  obscured  when  the  numerator  in  the 
abscissa  is  reduced  by  an  amount  equal  to  the  submerged  weight  of  the 
volume  of  sediment  displaced  by  the  embedded  object. 

From  these  data,  the  equation  for  computing  breakout  has  been 
determined  to  be: 


F  -  0.20qdAm,Ke-°'°0,M01,-2eo)  (7! 

where  Am<„  *  horizontal  projection  of  the  maximum  contact  area 
por  t  ■  0,  Equation  7  reduces  to 

F  -  0.81  qd  (8) 

For  t «  260  minutes,  the  breakout  force  requirement  becomes 

F  -  0.20qdAmi„  (9) 

It  Is  to  be  emphasized  that  relationships  such  as  appour  in  Figures  18a, 
ifib,  and  18c,  of  which  Equation  7  is  typical,  ora  based  on  only  one  type  of 
sedimont  (that  Is,  that  found  In  San  Francisco  Bay)  and  one  size  of  teat 
objects,  all  of  which  were  similar  both  from  characteristic  length  and  bearing 
load. 


THEORETICAL  ANALYSIS 

To  Introduce  the  theoretical  procedure  employed  in  this  study  It 
seen, 3  useful  to  pause  and  ask  why  and  In  what  way  Is  It  worthwhile  to 
proceed  with  the  development  of  complicated  models  and  procedures,  That 
question  is  best  answered  by  the  following  quotation  from  Whitman  (1934); 

"Highly  complex  theories  should  seldom  be  used  In  soil  engineering 
practice,  The  advent  of  modern  computers  does  not  change  this  conclusion. 

"There  Is  always  a  ba3lc  limitation  upon  any  computational 
procedure  that  Is  to  be  used  in  design  practice— the  user  must  be  able  to 
trace  a  clear  relation  between  each  assumption  and  the  result  of  making  this 
assumption,  In  the  case  of  soil  englnaering  dosign  there  la  an  additional  fact: 
It.  Is  almost  impossible  to  ascertain  the  actual  pattern  of  nonhomogonelty 
mat  exists  In  a  particular  soli  deposit, 
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"Sophisticated  models  and  procedures  are  useful  only  in  applied 
research.  The  end  objectives  of  any  development  of  such  models  and 
procedures  3ro  the  clarification  of  certain  puzzling  aspects  of  the  behavior  of 
soil  masses  to  bring  greater  unity  into  empirical  observations,  and  the 
determination  of  the  reasonableness  and  accuracy  of  simplified  methods  to 
aid  the  development  of  such  simplified  methods  for  use  in  practical  design 
siU'H'ions," 

The  essential  nature  of  the  breakout  problem  differs  from  the 
ordinary  footing  problem  (that  is,  prediction  of  the  ultimate  bearing  capacity) 
in  at  least  two  important  aspects. 

First,  the  ordinary  footing  oroblem  is  generally  analyzed  statically  by 
means  of  equations  or  charts  developed  by  Torzaghl,  Time-dependent  effects 
are  not  considered,  whereas  for  some  preliminary  breakout  tests  such  effects 
were  observed  to  be  Important.  Second,  the  loading  pattern  applied  to  the 
soil  muss  by  footings  Is  relatively  simple  as  compared  to  that  which  might 
exist  before  and  after  extraction  of  a  full-scale  submarine  or  deep-submergence 
vehicle.  This  Is  further  complicated  by  the  remolding  of  the  soil  In  the  vicinity 
of  the  object,  whereas  for  footings  tho  soil  Is  assumed  to  be  essentially  undis¬ 
turbed  The  degree  of  remolding  as  a  function  of  distance  Is  unknown, 
Moreover,  laboratory  techniques  for  determining  the  stress-r, train  relationships 
for  soils  under  compressive  loadings  are  well  established,  These  relationships 
are  useful  In  soil  engineering  practice,  since  the  tests  from  which  they  are 
derived  correspond  roughly  to  the  loading  patterns  Induced  by  footings,  On 
the  other  hand,  little  comparable  Information  has  been  developed  for  soils 
under  tonsil”  loadings,  Thus,  tho  symmetry  or  asymmetry  of  the  stress-strain 
relationships  for  soils  Is  not  well  established. 

In  view  of  these  differences  the  use  of  u  sophisticated  modol  for 
analysis  of  tho  breakout  process  seems  at  toast  partially  justified,  Tho  word 
partially  Is  sued  since  one  might  well  orgufl  that  It  is  difficult  If  not  Impossible 
to  model  mathematically  the  very  complicated  two-phase  material  doscr Ibnd 
generally  as  soil.  If  this  premise  were  accepted,  then  predictions  of  brookout 
behavior  would  be  required  to  rely  exclusively  on  experiment  followed  by 
the  use  of  relatively  simple  models. 

Thn  theoretical  procedure  used  in  tins  study  Is  based  on  a  systematic 
numerical  procedure,  developed  by  Harper  and  Ang  (1963),  for  determining 
tho  displacements,  atrolna  and  stresses  within  a  plane  continuum  wherein 
certain  regions  have  buon  strained  beyond  un  elastic  yield  limit,  The  material 
of  tho  continuum  Is  considered  to  bo  Isotropic,  olnstlc—por furtly  plain. c,  and 
thu  problomsaro  solved  for  continuously  Increasing  external  loads. 

The  numorlcal  procedure  Is  asaomlally  a  relaxation  technique  applied 
to  a  discrete  physical  model  composed  of  suitably  arranged  stress  pointu  and 
mass  points,  Once  the  external Iv  applied  loads  have  been  raised  to  a 
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sutticientiy  n ign  iovci,  the  mum  iiiyiiiy  ai, ua.ntivl  stress  point:;  of  th;:  mode! 
bogin  to  yield  and  flow  plastically.  The  initiation  of  yielding  is  determined 
by  the  Mises-l  It-ncky  yield  criterion.  Thereafter.  ylnMuri  me  in  ns  nrn  assumed 
to  obey  the  plastic  stross-strain  relations  postulated  by  the  Prnndtl  Ruoss 
theory.  The  numerical  method  is  presented  in  complete  detail  in  Appendixes 
C  through  G.  Only  a  brief  resume  Is  presented  In  the  following  paragraph, 

The  lumped  parameter  model,  shown  In  Figure  C-1,  Is  used  In  the 
procedure  developed  by  Harper  end  AnQ,  The  model  consists  of  moss  points, 
at  which  the  mass  of  the  material  Is  assumed  to  be  concentrated,  and  stress 
points,  which  connect  the  neighboring  mass  points.  Displacements  In  the 
continuum  are  defined  nniy  at  the  muss  points,  while  stresses  end  strains  are 
defined  only  at  the  stress  points,  T  he  stresses  at  each  stress  point  uio  rululod 
to  the  movements  of  the  four  surrounding  mass  points  by  using  a  finite- 
difference  version  of  tho  stress-strain  lews  for  the  continuum.  Boundary 
conditions  are  given  In  terms  of  either  external  loads  acting  on  tho  mass  points 
or  specified  displacamonto  of  those  mess  points, 

Figures  19  through  21  are  results  of  on  example  obtained  by  application 
of  the  numerical  method,  In  this  example,  It  Is  assumed  Ihnt  e  4  *4  x  lb-foot 
parallelepiped  has  ponotratad  a  distance  of  6  Inches  Into  n  soft  marine  sediment 
having  a  yield  strength  In  simple  tension  of  0,6  psl,  a  Poisson's  ratio  of  0,4,  end 
an  effective  modulus  of  elasticity  of  74  psl, 

Figure  19  depicts  the  propagation  of  plauilr.  straining  witn  Increasing 
force  levels,  Figure  20  shows  the  distribution  of  yield  rutin  porcontlloE.  for  a 
given  level  of  loading,  Tho  elastic  displacements  at  first  yielding  are  shown  In 
Figure  21,  Results  given  In  Figures  19  through  21  era  based  on  an  nrbltrurv 
loading  of  tho  boundary  men  points,  whirTi  duo*  not  nolle  coincide  with  the 
loading  pattern  In  the  field  tests,  The  application  of  uniform  loads  to  the 
mass  points  (that  Is,  land  Intensity  Is  constant)  results  h  unequal  movement 
of  those  mass  points,  which  corresponds  to  what  happens  with  u  flexible 
footing,  But  submarine  Keels  and  other  appendages «r« very  rigid  and  con 
respond  to  a  rigid  footing,  Thus,  the  results  given  In  Figures  19  through  HI 
must  be  considered  with  this  reservation,  Certain  difficulties  In  adjusting  the 
load  Increments  within  the  plastic  range  have  thus  for  prevented  the  processing 
of  dlsphcement-controlled  loading  pattern  problems,  l>>  summary,  no 
anomalies  ere  present  in  the  theoretical  approach,  even  though  only  the  gross 
feature*  of  the  breakout  process  have  thus  far  been  .diid'ed, 

In  order  to  substantiate  the  results  obtained  by  application  of  the 
numerical  technique,  a  photoolastlc  >  Uidy  id  the  stress-strain  pet  terns 
Induced  by  tunslle  loadings  on  various  boundary  geometries  wits  completed, 

The  details  ere  given  In  Appendix  IT,  A  umthnne  m hirer  having  n  modulus 
of  elasticity  of  600  psl  and  a  Poisson's  redo  of  0,45  was  employed  for  these 
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lusts.  The  tails  sunport  Iho  results  of  thu  numurictii  technique  in  ovory 

[]j  jnfiji1  w!ttii<<  tl'ii  elastic  whrrp  tlir  rt'*uj!t,i  nrr  c'tvnpnrat'!*-)  *1 '^n 
condition  of  plnsticlty  c.iinnot  be  simulated  by  the  photooiastic  method.  Thus, 
the  niimericnl  technique,  which  Is  applicable  to  both  the  elastic  and  plastic 
regimes,  Is  much  moru  versatile  than  tin  photoelostic  method  and  can  bo 
applied  to  a  wide  tango  of  problems  not  necessarily  llmltod  by  laboratory 
equipment  or  Iho  construction  of  apodal  modal*.  Tim  prostint  vorslon  of  thu 
immnrioal  method.  hnwavor,  Is  limited  to  a  two-dimensional  trontmnnt  ol 
iho  problem. 


COMPARISON  BETWEEN  THEORETICAL  AND  EMPIRICAL 
PROCEDURES 

An  example  Is  selected  to  compare  the  lasults  nf  the  theoretical 
procedure  with  (hose  of  the  empirical  formula.  Data  from  the  photoalastlc 
modal  arn  also  used  to  estimate  the  breakout  force, 

Let  us  assume  that  the  keel  of  the  oudeer  research  submarine  NFM  Is 
burled  6  Inches  In  a  marine  sediment  similar  to  that  In  Sen  Pranclyot)  Bay, 

This  soil  Is  saturated  day  and  can  bo  assumed  lb  have  a  Poisson's  ratio  of 
0,40  and  a  yield  shear  strength  u I  0,0  pul,  Let  us  further  assume  that  the 
effective  modulus  of  elasticity  of  the  soil  is  74  pal,  which  Is  etpial  to  the 
average  modulus  of  elasticity  of  the  tested  soil,  This  value  Is  obtained  through 
trlaxlal  tests,  The  keel  is  considered  as  a  parallelepiped  Oh  feel  limp  and  4  feet 
wide, 

The  breekout  force  requlretf  fur  an  Immediate  pullout  Is  first  iialoulatad 
by  the  analytical  method  given  In  detail  In  Appendixes  C  through  Cl .  The  soil 
properties  end  the  nhjuot  geometry  are  mitered  as  Inputs  to  the  computer 
program  according  to  the  format,  The  force  Inclement  Is  then  lelcuted  to  give 
the  desired  accuracy,  The  final  force  which  t  mixes  dm  soil  in  yield  anti  the 
nh|ei,(  to  totally  detach  (tom  the  bottom  Is  the  hieaknut  force,  This  fume 
has  been  found  hum  Ihecomputai  output  to  he  slightly  laiger  then  dl.O(X) 
pounds, 

If  the  soli  is  considered  e  pule  elartlc  material,  then  the  baeknul 
force  may  be  estimated  Iron  m  pholoe'asllu  model  dale,  based  on  model 
sillily  theorem,  we  have 


where:  m  »  modul 

F  -  breakout  form 

L  “  lunglh  ot  purellolepiped 

B  “  width  of  poralioloplpud 

Tm*«  "  maximum  shw.tr  stress  or  strength 

In  this  case 


"  6, 14  pounds 

L 

*  9ft  feat 

«  0,7(5  Inch 

B 

■  4  font 

®rn 

■  0,26  Inch 

“  0,6  p&l 

)m“  11,6  pal 

Substituting  these  values  In  the  oquutlon  we  obtain1 


F  -  J,14(95).(4)(144K0.6) 
0,76(0,261(11,8) 


03,600  pound* 


The  Poisson's  ratio  for  the  modal  material  Is  0,46,  The  photoelastlc  model 
study  Is  discussed  In  Arvundl*  H, 

The  brookout  force  Is  also  cttlculatml  (runt  the  empirical  formula 


F  -  0,20qtlAmMo‘H,','«1 


whom  qu  ■  2,86  [1  +  (B/L)]qu 
In  this  problem 

AmM  -  1)6(4) (144)  In.!*  “  64,1)00  ln> 
t0  -  260  minutes 

t  “  0  min  (Inst, in,  breakout) 

R  0,0664 


0-4  (net  (width  of  paralleloplpod) 

L  -  Ob  (out  (lonuth  of  parallelepiped) 
qM  "  6,6  psl 


where  qu  Is  the  uni  onflned  compressive  strungth  raid  Is  usiumod  to  be  equal 
to  the  yield  strength  In  simple  shorn. 
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Substituting  these  values  in  tne  empirical  iummloc,  wg  have 


qd  -  2.85(1  i  (4/95)]  (0.6)  -  1.79  osl 

Thus  F  »  0.20(  1.79) (54,900)e"°'oos4<0 ' 2601  -  80,000 pounds 

Thus,  comparable  predictions  of  the  breakout  force  for  zero  elapsed 
time  according  to  the  theoretical  and  empirical  approaches  are  91 ,000  and 
80,000  pounds,  respectively. 

The  theoretical  procedure  considers  neither  the  effects  of  remolding 
or  the  increase  in  strongth  due  to  consolidation.  In  the  San  Francisco  Bay 

tests  (upon  which  the  empirical  constants  were  determined),  the  maximum 
bearing  loads  were  very  high,  much  higher  than  the  capacity  of  the  soli  near 
the  surface  to  support  3uch  loads,  Thus,  the  soil  in  the  Immediate  vicinity  of 
the  test  object  wus  disturbed  and  remolded  to  a  considerable  degree,  More- 
ovai,  the  object  penetrated  a  certain  distance  until  the  bearing  loads  wore 
reduced  to  a  level  within  tho  capability  of  the  soil  to  support  the  imposed 
loads,  There  i»  a  natural  Increase  In  strength  with  depth  due  primarily  to  an 
Increase  In  bulk  density  and  a  decrease  In  watar  content,  However,  the 
strength  of  tho  soil  Is  also  affected  by  the  presence  of  the  object  in  two 
opposing  ways,  One,  referred  to  previously,  Is  the  reduction  In  strength  due 
to  remolding,  The  other  Is  the  gain  In  strength  due  to  consolidation.  Both 
effects  occur  on  vastly  different  time  sculos.  The  loss  In  strength  due  to 
remolding  takes  place  Instantaneously,  whereas  tho  gain  In  strength  due  to 
consolidation  Is  a  long-term  process  depending  Initially  to  a  large  extent  on 
tho  permeability  of  tho  soil, 

It  seems  likely  that  for  a  given  object  two  worst  situations  are 
possible,  One,  termed  tho  shallow-penetration  case,  occurs  when  tho  soli  has 
a  high  shear  strength  which  Is  almost  but  not  quite  matched  by  the  Imposed 
bearing  loads,  This  ensures  ti  dose  bonding  of  the  object  skin  surfuce  to  tho 
sedimentary  layer  without  Inducing  a  strength  reduction  In  tho  soil,  The 
other,  termed  tho  (loop-penetration  case,  occur:,  whun  poriutrotlon  hus  been 
so  deep  that  tho  volume  of  dlaplucud  soil  becomes  sufficiently  large  so  tin  to 
completely  dominate  the  breakout  process.  The  latter  situation  Is  not  to  bo 
confused  with  tho  volume  of  motorlal  lying  be  tween  the  fulluro  surface  mid 
the  object  boundmy,  which  Isa  function  of  the  gross  dimensions  of  tho 
object,  We  are  concerned  here  with  a  given  geometry. 

In  the  theoretical  approach,  tho  computational  scheme  permits  loads 
of  any  magnitude  (o  be  applied  to  any  or  all  m«SB  points,  Again,  In  tho  cited 
example,  equal  loads  wore  applied  to  the  mass  points  locutod  on  the  boundary 
goomatry.  Thus,  unequal  relative  movement  Ixi’woan  the  mass  points  on  the 
boundary  are  permitted,  whereas  In  fact,  such  unequal  movements  are  realized 
only  for  flexible  membranes, 
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FINDINGS 


1.  A  numerical  method  of  predicting  strains,  stresses,  and  displacements  in 
an  elastic — perfectly  plastic  medium  subject  to  loads  applied  to  an  arbitrary 
boundary  geometry  was  found  to  be  useful  in  developing  a  theoretical 
prediction  of  breakout  forces, 

2.  A  complicated  computer  program,  which  uses  a  lumped  parameter  model 
of  the  material  and  an  iterative  technique  to  obtain  solutions,  was  found  to 
bo  an  integral  part  of  the  theoretical  procedure.  Tho  program  requires  use  of 
a  high-speed  large-memory-capacity  computer, 

3,  Tho  computational  procoduro  traces  tho  development  of  the  stress  and 
displacement  fields  in  on  olustic  -porfoctly  plastic  material  under  conditions  of 

plane  strain,  with  specified  boundary  conditions  and  force-controlled  loading, 

4,  Results  yielded  by  the  computational  procedure  were  found  to  be  verified 
by  separate  photoelastic  studies,  at  least  within  the  elastic  range, 

5.  Data  from  breakout  tests  with  large  specimens  in  San  Francisco  Bay  were 
found  to  develop  tho  following  empirical  formula: 

F  -  0,20  Ami„  qd  t*000#40*1  •  J8°l 

The  geometry  of  the  breakout  object  seemed  to  have  relatively  little  effect  on 
the  breakout  force, 

6,  In  a  portlculur  example,  the  breakout  force  requirement  was  estimated  by 
the  theoretical  procedure  to  bo  91,000  pounds  and  by  tho  empirical  procoduro 
to  bo  about  B0, 000  pounds, 


CONCLUSIONS 

1,  The  ocean  bottom  breakout  force  of  on  object  of  simple  geometry  cun  be 
ustlmutod  by  moons  of  an  analytical  method  thut  uses  numerical  calculation 
by  high-speed  computers,  Tho  method  tokos  Into  account  tho  plastic  behavior 
of  soli  beyond  the  elastic  strain  range, 

2,  Tho  following  empirical  formula  may  bo  used  to  describe  tho  breakout 
force  for  an  ocean  bottom  soli; 

F  *  QA1„Hqt,rM’-'o| 
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recommendations 


1,  Whenever  possible,  engineers  should  use  the  analytical  method  outlined  in 
this  report  and  the  appendixes  to  estimate  the  breakout  force  before  any 
actual  salvage  operation, 

2,  When  a  computer  is  not  available,  the  empirical  formula  should  be  used  to 
determine  the  breakout  force,  The  constants  Q,  R,  end  t0  can  be  derived  from 
a  limited  number  of  in-situ  field  test  data, 

3.  More  field  data  should  be  collected  to  verify  the  analytical  method  for 
various  object  geometries  and  soil  sediments, 

4.  More  research  should  be  conducted  Into  methods  for  reducing  the 
breakout  force, 


NOMENCLATURE 


9 

A  Horizontal  projection  of  the  contact  area,  in. 

Am(ll  Horizontal  projection  of  the  maximum  contact 
area, in,2 

B  Width  ol  breakout  object,  foot 

b  Numerical  constant  used  to  force  C  “  C  for  very 

large  t 

C  Collusion,  or  nltnrnutivulv  a  moasuro  of  tho  vans 
shuor  strength,  psl 

C  G ffoctivo  average  cohoslon  along  tho  falluro  surface 

at  tho  Instant  of  breakout,  psl 

F  Broakout  forco,  pounds 

k  Constant  which  Is  o  function  of  object  size,  object 
tthnptt,  time  duration  ol  applied  force,  rate  force  Is 
applltid,  soil  sensitivity,  and  tho  ulopaod  time  which 
the  object  has  buun  In  place  after  the  Initial 
disturbance 

L  Length  of  breakout  object,  feet 

Q  Constant 

11,1  Average  supporting  pressure  provided  by  the  sell  to 
maintain  the  embedded  ob|oct  In  static  wcnilllb- 
rlum,  psl 

linconf'nml  compressive  strength,  i>sl 

Ft  Slope  ol  thu  "failure  lino"  when  log  IP/gA)  Is 
plotted  versus  time,  t 

S  Dngnm  ill  soil  Minsitlvitv 

t  I  line  allowed  fur  breakout,  minutes 

t,  Miifiiruni  e  limn  minimi  In  thu  thixotropic  behavior 

if  a  mnieilal  In  mgalnlng  e  Muled  pen  Milage  of  its 
iliiuigtli  nftie  initial  dlsturlMiu  n,  minutes 

t(1  >ti<f(ii in ii  it  time,  mlnutiiii 

T n  i .An •  111111111  slums  xli i'h‘i  i u  strength 
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Appendix  A 


SITE  SELECTION  AND  LABORATORY  TESTS 
by  M.  C,  Hlronaka 


Tho  field  tests  worn  carried  nut  In  San  Francisco  Bay  about  BOO  yards 
southeast  of  Pier  3  at  Hunters  Point  (Figure  A- 11.  Sen  Fn incisco  Bey  was 
selected  because  it  Is  the  only  location  on  the  West  Coast  In  relatively  shallow 
water  where  there  exist  large  undisturbed  deposits  ol  line  grain  sediments 
whose  physical  characteristics  nro  reasonably  similar  to  those  of  the  deep 
ocean  basins, 

Tho  specific  site  tanno  which  Is  relatively  free  of  the  shells  and  debt  Is 
which  are  generally  found  adjacent  to  piers,  wharves,  and  sea  wells,  and 
which  make  such  locations  unsuitable  for  field  tests,  Figure  A*2  shows  the 
depth  contours  at  the  site  and  the  soil-profile  locations  along  which  soil 
cores  wore  retrieved,  Figures  A>3  through  A-B  are  schematic  elevations  of 
the  soil  layering  based  on  a  detailed  examination  of  the  soil  cores,  Table 
A-1  gives  statistical  Information  on  the  laboratory  soil  tests  for  certain 
selected  cores, 

Finally,  Figure  A-B  presents  a  typical  sample  of  the  results  of  tilaxlal 
tests  of  certain  soil  samples,  Voices  of  cohesion  obtained  from  the  biaxial 
tests  Hppeerotl  to  agree  with  corresponding  values  obtained  *rom  vane  shear 
tests,  The  bulging  of  the  soil  sample,  which  occurred  In  the  latter  stage  of 
the  trlexlal  tests,  tends  to  make  tho  remainder  of  tho  biaxial  test  results 
somewhat  less  than  meaningful, 


33 


I  iihli  A  I ,  i  ,-sl  I  lest  ills  I  mm  Sc  hi  led  (‘me  Sum  i|  ill's 


* 


Cut#  Sumpi*  JO 

(l\:v  *i\gln  ■  0?.tn,  radial  illaJnm «  «  tI0 1  Ml 


I'f  (lixtrtlui 

I'lllllvlll  III!,) 

l?.1tt 

3437 

1 

40.01 

73./0 

114  87 

0(1-00 

I0HII1 

OiiiV  w«t  Uithftlty  <1**0 

01 

Oft 

08 

UK 

UK 

04 

80 

01 

Vh»>  il  Mil  ttlcnutli  li'*M 

(V) 

11? 

104 

1 '2 

101 

108 

INI 

01 

MumoliKn;  miiihgtii  (|>*I I 

11 

13 

3.1 

18 

34 

?:l 

34 

III 

Rnnilllvll" 

M 

U.1 

4,0 

n.i 

n,:i 

7,:i 

(>.7 

:i,s 

OHill'idl  waluf  rniMiml  W) 

lot!  ;i 

u:».  i 

111,4 

Uli.o 

110,3 

8/0 

711  4 

7.1,11 

BtuttOlo  uilvlly  nl  lolldl 

a.ui 

3,03 

3,0! 

2,00 

.7.07 

2,00 

?,04 

2,03 

Dry  (IhmiIIv  IirII 

44 

SO 

M 

40 

03 

47 

04 

02 

Vulil  mile 

2,  U 

3,10 

3,02 

3,32 

3,32 

2,40 

2,00 

2,11 

ftoMUly  1%) 

7.1,3 

00,7 

00,0 

EH 

(11,0 

71,3 

07,3 

Mtumivd  void  min 

3,11.1 

3,44 

3,1  ! 

2,40 

2,37 

2,b7 

1,00 

1,04 

Liquid  llmll  (HI 

7(1 

04 

02 

7il 

H3 

04 

04 

I’lUIKl  llmll  IHI 

40 

l‘S 

ft# 

30 

HI 

,10 

ClMlIuily  IimIm 

13 

31 

,14 

43 

.14 

34 

Llnuldlly  Indo* 

307 

1U0 

Iftft 

101 

too 

1,14 

1.14 

1 38 

U.imiwvitlim  HHlMN 

A3 

,40 

■iM*4M4'  IUkJ 

,4? 

44  W  A  OM, 

,07 

,N) 

.1)0 

*4  h_a**L_ 

,40 

,40 

ttmilimpli  M 

(Mil*  mgl*  ■  0fl,0l,i  fttllil  illltlitt*  «  741  III 


•HutMMlM 

— * 

IliIMVl!  Ill*,) 

12,10 

34.27 

30.30 

0003 

72,1ft 

04.07 

00  DO 

1  Oft, HI 

Hulk  m*i  ikmiliv  Hwll 

04 

02 

07 

Oft 

03 

Oft 

um 

on 

V«n»  iiiott  iinmglii  (noil 

Ml 

02 

10.1 

Iftl 

141 

II? 

407 

HI 

llttmuldwl  (imiylli  lidill 

11 

17 

1ft 

IH 

24 

20 

411 

II 

IhmiltMly 

ft, 3 

4,0 

ft,  7 

0,3 

U,0 

,1.0 

0,0 

4, II 

(O,0l84l  mill  (1881*81  IHI 

100,0 

1UU.li 

IU4, 7 

03,7 

OH, ft 

011,1 

(13,7 

(ll.ft 

HiMH'IIIV  0I4VIIV  Hi  kllllH 

2.  Ml 

3,00 

2,71 

2,73 

2, DU 

2,00 

Oittft 

2,72 

(liy  iinfiikllv  IimII 

40 

4ft 

LI 

HI 

47 

0? 

111 

1)4 

Vulil  0*118 

2,07 

2,01 

2,30 

2,10 

2,03 

2,2ft 

1,71 

2,1? 

Cuoiilly  IHI 

12,1 

72,3 

DU,  7 

no.  7 

71,7 

00,2 

0.1,1 

MU) 

0418111*1  yulil  mild 

2,0'l 

2  07 

2,40 

2.70 

2,  Oft 

2,31 

I,IW 

2,10 

Utiulil  limit  IH) 

00 

72 

H7 

1)11 

77 

04 

Ml 

Oft 

ftlMlI,,  llmll  IHI 

.1.1 

It 

W 

,11 

;t« 

20 

:ut 

:th 

CImIIiHv  Imhu 

'W 

:to 

30 

14 

:w 

?ft 

37 

3,1 

Mqillllllv  Kill** 

2(01 

170 

172 

Ift4 

140 

222 

Hft 

1,18 

L'(mi|il**lii,i  Uni** 

.07 

,11ft 

,1)2 

,0.1 

,W) 

,1? 

.  .  •  S'. 

,11? 

CmmIIiiiiimI 


M 


t  lllllr  A  I  OlhtilHI'Ml 


If 

IISiMi  mn|li>  ■  III)./1,  t  rtillnl  tlltlitttt  it  “  -11)11111 


I'll  lIVKl  II-  1 

MB 

0:1 

U  Hi 

>4  Zf 

in  hi 

40-01 

non:: 

//  ft) 

N4.it/ 

HulV  Alii  ilnitiltv 

H4 

»' 

114 

IHl 

no 

in 

lift 

M 

Vtnit  tlmw  «l‘i>nglh  i|»ll 

tt:i 

in 

iM 

11/ 

i  i.i 

»4 

1W 

Ml 

J  \»»tiu  iIiIik  1  til im(| t It 

\n 

in 

14 

1 ! 

in 

14 

Zl 

i;i 

IftMUlllvlly 

XU 

•\M 

n  A 

U.W 

l,\ 

ft. <1 

U.t\ 

,i(i 

OtlgliHII  v.Mxl  mlittlll  (Hi 

«/  ,i 

104,0 

01, 1 

no, ii 

114,0 

00,0 

Oi.tt 

00,0 

HimllH  guvlly  til  tuiltlt 

II, iM 

4,0/ 

/,11ft 

/,0ft 

z,ow 

0,0/ 

11.00 

Z.Pfl 

(Hy  ilki'tiiy  (fvrll 

no 

4ft 

40 

4/ 

4/ 

411 

0/ 

DO 

Villtl  Mill' 

l,W4 

/,/  4 

Z„(Z 

it, 04 

1,04 

1U)0 

11,1/ 

It,  14 

Otllllllly  (HI 

00,1) 

/II, II 

Z0,1 

T4,n 

/ft,/ 

/)/,/ 

00,0 

70,0 

tOMntiliKl  viiltl  ralln 

UN 

/,/ft 

0,4)1 

,1,01 

,'1,00 

It, 00 

it, in 

it, 04 

l,H|Hltl  limit  (HI 

40 

/I 

«/ 

ON 

0/ 

VI 

0* 

ZU 

HlMIlii  limit  IHl 

IN) 

11/ 

11 

41 

41 

Ml 

It! 

,10 

I’Ulllully  lllilu* 

III 

DO 

14 

4/ 

40 

4/ 

40 

41 

Lltliiltllly  itiilt* 

lift) 

IN/ 

lift) 

10)1 

100 

1.10 

140 

U» 

Cninoiiikilmi  iihIm 

,11 

,li|) 

,01! 

/a.-  x  a.a 

zi 

,Z0 

i>8«  '^1 

,00 

,01 

,nst 

U'liitiimi*  tv 

IIMItl  tligfti  “  Util",  I4tlltl  llltlMMVM  «  41/11) 


k— i-.,VA*tU»— .  «.**■»_ :-.i;  (in  i:  ir.w 

tftft|tHflfti« 

i  — - 

IiMkI 

v«!  lift,) 

.is  .  : 

l/illi 

444/ 

14), W 

40  01 

14)0,1 

z/Zh 

tM-Nf 

BOftU 

Hull,  *y|  lUntlly  Mu  ll 

11/ 

ftl 

«/ 

ftl 

Ml 

ftl 

111 

Ult 

Vtnt  iftn«  iliMiiyin  Mull 

04 

/O 

HO 

04 

lllft 

144 

100 

IUI 

llmimlilml  ilivngili  l|i*ll 

ft 

1/ 

'4 

tit 

,  // 

'/II 

IN 

11/ 

ftWitllWIIV 

0,0 

0,0 

0/ 

0,4 

4, ft 

M 

11,0 

Ift, 4 

mill  inti  wtwi  tiiiitiviti  IHl 

l.l/.N 

104,0 

10, Ml 

ftl, II 

141,11 

0:1,0 

«4,0 

ft/,' 

Nimt'ilin  tiruvlly  nl  tnluit 

it  111 

/,I1> 

/,0li 

/.till 

4, ON 

/.lift 

/,OII 

/,zo 

Illy  llMII'lli  1|H't) 

,MI 

40 

40 

411 

411 

ftl 

4/ 

40 

Vnltl  itlln 

.till 

/,ft/ 

/,/'/ 

/,4/ 

4,44 

J,4ft 

/hi 

11,00 

1*ill,*lllf  IHl 

//,ll 

i/,(i 

/l, I 

10  / 

Mil 

141,1 

t/,o 

/1, 4 

lUlilltllHl  Vtllll  Itlln 

M'  1 

/,/  * 

/no 

Mil 

/.iftl 

/  /I 

/.fti 

Z,4U 

1  I'lnlil  ftii'll  IHl 

UO 

/ft 

it/ 

// 

til 

nu 

HI 

/ft 

I'lttlli  lln,i|  IHl 

At» 

Ih 

HI 

in 

41 

14 

40 

HI 

Iftttlli'lly  In, Iti 

INI 

40 

40 

m 

1/ 

10 

4:1 

40 

1  ImuIiIM V  null1* 

1110 

1// 

inn 

U4I 

10ft 

140 

l/ft 

III 

Cn,ii|iiu»iiMi  imii,  > 

.1/ 

■ 

,llU 

,0ft 

,0ft 

.04 

,01 

,00 

,0/ 

.11) 


Ciiiiliftiiytl 


InliloA  t.  Contliuiml 


Cot*  Bumpli  tw 

|l’(il*r  miglr  »  UP,Ou;  rsdlll  iHitnnco  ■  602  III 


I'roimttlHt 

liltnii'ill  tin.) 

•;4.;‘7 

14)  OH 

Tin 

100.111 

Hulk  Wl»l  linniUy  (|kM 

HV 

Kill 

in 

112 

HI 

M.1 

00 

»4 

Vmiu  tli ««  iLwtgll,  (pill 

U.l 

Ub 

70 

100 

124 

12U 

1 24 

131) 

Homolttod  llfunglh  ||i»t) 

»:i 

III 

16 

2U 

20 

31 

46 

10 

S*n»Hlvlly 

ii.d 

0,1 

4,4 

3,7 

4,3 

4,1 

2,7 

13,6 

Orlglnll  W1O01  mmliml  IK) 

1 011,0 

bU,1 

107,6 

64,6 

04,7 

67,6 

62  3 

83,6 

RlHtulllu  grtvlly  id  tulltli 

2,08 

2,00 

2,66 

2,66 

2,66 

2,67 

2,66 

2,66 

Dry  dmulty  (poll 

40 

67 

44 

60 

47 

60 

63 

61 

Void  min 

a,  oa 

1,02 

2,70 

2,34 

2,63 

2,36 

2,1b 

2,24 

Porotilv  (X! 

711, 4 

60,3 

73,6 

70,1 

71,7 

70  3 

66,2 

66,1 

B«UirylMl  void  rtllo 

a,  70 

i.btt 

2,67 

2,26 

2,62 

2,34 

2.10 

2.21 

Liquid  limit  :WI 

74 

43 

63 

60 

112 

77 

72 

76 

Pldllu  limit  (HI 

Ml 

22 

40 

33 

36 

36 

36 

16 

Ptantoltv  IwIm 

30 

20 

41 

36 

44 

42 

30 

'49 

Liquidity  IruM 

176 

IU0 

167 

143 

126 

124 

127 

123 

COh  i|'f»WlOII  ImlMK 

,t17 

,20 

,66 

,62 

,66 

.61 

,66 

,66 

Ml) 


I— JU-  I, - I 

0  100  200  300  HorltonOI 

0  “to  JO  30  Vlfllcil 
Heal*  (III 


toll  oliy-* 
mitllum  dim  dliy 

dim  ttlay 


Plflur*  A-3,  ftohwiiMle  of  loll  liyirliifl  *lonu  pruflli  A, 


SFH  ilitll  li'agmtnli  ind  Mnd 
PT  ally  nut  wnu  ihill  digmmli 
HJ7  id  ill  Irigminli 
(U  ai«y 


iKinom  ui  aoii 


mini  loy*r  luw  wont 


0  106  ilk)  SHK)  HwImmiI 

(i  |6  sb  "Jo  VlMIlil 

lull  111) 


mill  li»|m»ini»ml  wml 
illV  imltumi  ihill  Irigiiiinli 
mill  Iriiminti 

LLl  illy 


litiHiiitt  ul  Kill 


Pl(|un  A'4,  Boh  wiinf  III  of  mil  mynfl««  ulmm  imifllo  If, 


il  f  So . iftrt  siVi 

Wl>  Mltti.  «i  tel  I  4 

0  ii)  ilo  ,Vi 

Nlllll  (Ill 


Him  yliv 
iimlliim  linn  I'liy 


llillKlllllil 

Vii  Hall 


I 


wull  Hi|iiiiiiit  iiiil  mill 
,'liy  nul  until  mill  Hlyminlt 
IAI  Pull  Hmniinli 

El  »i»v 


Him  illy 


lllllHlIIII  dim  yliy 

'O'  Iwiiiim  itf  mil 


hyiiri  A  D  Huliimillu  ot  Mill  liyirlny  ilonu  ihiiIIIm  (!. 


Append*  B 


LITERATURE  SURVEY  OF  BODY  SHAPES  AND  DIMENSIONS 


h-  r>.  J,  Muga 


A  thorough  rnvlew  of  tha  llteruture  on  existing  body  shapes  and 
climnnnlc-nn  of  possible  Inteiaat  to  the  breakout  foroo  problem  was  conducted 
for  thii  purpose  of  establishing  specific  guidelines  lor  designing  the  experi¬ 
mental  lest  piogrmn,  One  large  dess  nl  hn  ly  shapes  nemely  anrhms,  was 
excluded  from  this  study,  It  was  felt  that  the  Irregular  nonformall.etl  shapes 
presented  by  anchors  constituted  a  complete  study  In  Itself  and  that  there 
would  be  little  to  gain  from  the  attempted  extrapolation  of  data  available 
from  tests  on  anchor  shapes  to  othur  shapeii, 

Although  It  was  known  that  the  shapes  of  certain  classified  objects 
would  be  of  Interest,  they  also  wore  excluded  In  order  to  maintain  a  wide 
dissemination  and  distribution  ot  the  results  of  the  test  program,  The  study 
war.  focused  on  two  categories  of  body  shapes,  The  first  category  Included 
fu'hslra  U,  8,  Navy  submarines,  (he  various  types  nl  whluh  are  indicated  m 
T  side  II*  I ,  All  of  these  vessels  have  essentially  a  cylindrical  hull,  However, 
r,t  least  one  of  the  target  submarines  >881)  has  a  prismatic  ractangulai  keel 
attached  to  the  underside  of  the  hull,  The  second  category  Included  the 
gmwlng  family  of  deep  submergence  vehicles,  summail/ed  In  Table  ll-tf,  1 1 
Is  to  be  noted  that  toe  OOLI’HIN  (AC188*tMM)  Is  listed  in  both  catugntlas, 
However,  date  on  the  DOLPHIN  pertinent  to  thehteekout  force  problem 
did  not  eppaar  In  tne  open  published  llteratuie,  Of  the  othei  vessels,  the 
MORAY,  NILI,  and  .»K*lhseem  tu  fell  Into  the  leryei  dess  id  deep  sub* 
mergence  vehicles  having  appiuklmalely  the  same  contact  aina,  hot  some 
vessels,  such  as  the  AlllMINAUT  and  ALVIN,  data  on  stdface  contact  area 
was  not  leadlly  available, 

The  shapes  o*  the  bottom  contact  members  of  tlm  deep  submergence 
vehicles  vai  l»:l  from  ellipsoids  and  prisms  (rectanyulie  and  cylindrical)  having 
liiiyn  lengllvpndlamatni  t iitlntt  tu  morn  compact  bodies  such  es  cones  and 
Hphi'ies, 

I  loin  |bu  foittgulng  tribulations,  the  bash  shapes,  trims  and  welglds 
went  selected  for  |be  tost  pntuitim,  Nl/eit  and  weights  warn  suite  l«d  on  the 
bests  of  ID  beating  load  similar  to  that  which  might  be  experienced  in  the 
worst  possible  situation,  U'l  suitace  contact  aiee,  and  III)  available  weight* 
handling  capability, 

Dimensions  ami  additional  descriptions  nl  tlm  lest  objects  are  shown 
in  I  able  IK), 
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THtOntTICAL  ANALYSIS 
Ijy  0,  J.  Mugs  (iiul  W.  fj.  Atkin* 


INTRODUCTION 

I  hn  i  nsl  ul  lull  si  nlt>  bumt-unl  insls,  with  \  ;>fyiui  Inst  s|in<  inmns  mul 
si  til  |iiii|ii'i  Ill’s,  mnkns  II  idlin'  live  tn  lumnilnli'  o  mss  i*>|n>Msivi’  ii  ml  Inn  I  ul 
DH’illi  tli'||  I'l I'llki nil  1 1 n i  i’S.  I  his  (i|i|i"inll  >  1 1 ’In il  1  s  i 'u  n  Ihi'i ill'll'  ill  s| i nls' 

which  Is  hinted  on  n  tui'tHttlcAl  method  ol  iiMidltnlng  hot  him  littmktHil  lot'Win. 
It  Is  Ihtnmlml  Hint  eoidldtmon  urn  ho  dovelnjintl  in  thn  mothml  by  CHlIhrnllnu 
thn  (tnitlyllt ul  mnitol  tun  Milntivoly  snoill  numbm of  Itilloiciiln  fluid  lost*, 
t  hmooltm,  thn  rtnnlytUinl  method  mm  lm  timid  to  (mullet  liimikout  tindof  tl 
vmy  wltlo  ntngo  of  elrtiuiMtmum*, 

I  ho  following  sections  imistinl  In  dotml  iioillntmt  bmikgiimncl  lofijc* 
motion,  ilHsmibo  thn  numntkiol  inotliod,  develop  thn  oonMlhitlvo  onoollotis 
(mm  l ho  tliooiy  of  omtmilly  iilosllc  sulldn,  ond  liullcnfo  bow  Ihnsti  munitions 
con  bo  oomhlned  to  inovldo  n  svstnmotir.  mimmlmil  solution  of  uohtulnnd 
lihistk!  (low  omhloms,  (' (holly,  Homo  mutmplti  imihloms  mo  givnn  to  demon- 
sttntn  tlumseol  thn  oomontm  iimgimn, 


UACKQROUNt) 

I  lm | mi  mul  Ang  (  Hltllt)  Inwo  developed  o  syshmutllt:  niimmldttl 
piocmluie  Ini  tloloi mining  llm  dlsplm  ements,  *holni,  ond  siiesims  within  n 
pittite  cohtlhoum  wlmioln  mil  tnln  legions  loivo  boon  shinned  Imynnd  oil 
oloslln  yield  limit.  As  ndo|tbiil  lot  imu  In  this  study,  'ho  pincoiluin  Is  lestili'ted 
In  stotic  imiblonts  with  itlnno  sum  t  nmtllilnns,  I  ho  iniilotiol  ul  tbo  mmtliunim 
Is  nutkldorml  tn  lm  Isnhopie,  olosbc  iioiltmily  plnsth  mul  tlm  imihloms  mo 
solved  lot  coinimmiislv  Inc.iooslnu  osloliml  Inmls,  Unhnnllng  tium  o  plmtu  oil/ 
Mhlinmi  stoto  Is  mil  t  uiuldnind, 

I  lol|mi  oiii I  Ang  siolo  llod  "llm  nummli  ol  piu< eduin  Is  onswilinllv  o 
tnlilkilllnn  lot  Imh |un  oigillod  hi  o  t llm  mho  iihyslt  ol  mndol  1 1 n 1 1| k »um I  ul  milliitily 
itilungml  stioss  |  mints  ond  moss  |mlnts,"  ( tin  o  thn  ns  Inn  ml  holding  hits  I  mint 
inlsod  In  o  sntllclonlly  high  lovol,  tlm  Minn  highly  sliossint  inghois  ul  thn  no ••  ml 
Imgin  In  y  lull  I  mid  flow  tilosth  oily,  I  ho  Inltlotlun  ul  yielding  Is  iloloi  initmil  by 
llm  MIuiih-I  lonrky  yield  'ilhuiun,  I  hot  mil  lm,  yielded  loijiomi  nm  ossnmod  In 
obey  tlm  plmtllc  nIiosm  stmln  inloliuns  iiiulnlolnd  by  llm  I'lomltl  llmiss  tlmory, 


A! 


Use  of  the  discrete  model  permits  a  problem  m  continuum  mechanics 
with  an  infinite  number  of  degrees  of  freedom  to  be  replaced  by  a  correspond¬ 
ing  problem  in  particle  mechanics  having  a  finite  number  of  degrees  of  freedom. 
The  basic  odvantoge  of  such  a  technique  Is  that  it  makes  possible  the  solution 
of  problems  not  easily  solved  by  mathematical  analysis,  particularly  problems 
involving  partial  loadings  and  complicated  boundary  conditions.  This  Is  of 
considerable  Importance  to  tho  breakout  force  problem  since  many  different¬ 
shaped  objects  are  to  be  retrieved  under  a  variety  of  embedment  conditions. 

The  basic  disadvantage  of  the  use  of  a  discrete  model,  as  pointed  out  by  Harper 
and  Ang,  Is  its  very  finiteness;  that  Is,  stresses  and  displacements  are  determined 
only  at  a  finite  number  of  points.  He  frequently  the  finite  model  can 
furnish  only  «  rough  quantitative  measure  of  the  true  but  unknown  solution 
In  the  continuum. 

To  gain  somn  notion  of  the  accuracy  of  the  model,  Harper  and  Ang 
solved  a  problem  in  plane  elasticity  using  an  analytical  solution  and  compared 
the  results  to  those  from  the  numerical  solution  using  the  model.  The  results 
from  the  two  solutions  agreed  very  well. 

Whitman  (1964)  has  considered  the  use  of  the  numerical  method 
developed  by  Harper  and  Ang  for  analyzing  two-dimensional  and  three- 
dimensional  boundary  value  problems  Involving  soil.*  Whitman  concluded 
that  immediate  progress  with  the  problem  of  multldlmanslonei  contained 
plastic  flow  would  result  from  application  of  the  numerical  method,  and 
that  the  use  of  the  theory  for  a  perfectly  pustlc  solid  permits  a  first  step  In 
the  development  of  procedures  for  analyzing  contained  plastic  flow  In  soils. 

Christian  (1966)  has  used  the  numerical  method  developed  by  Harper 
and  Ang  to  Investigate  the  stress  and  displacement  fields  for  an  elastic  porous 
-  material  under  conditions  of  plane  strain  with  rectangular  boundaries  and 
force-  or  displacement-controlled  loading, 

Whitman  and  Hoeg  (1965)  have  examined  the  performance  and 
accuracy  of  the  Harper  and  Ang  numerical  method  and  used  It  to  analyse 
the  development  of  the  plastic  zone  beneath  a  strip  footing  resting  on  an 
elastic- perfectly  plastic  foundation  material,  The  computed  results  foriiie 
flexible  and  rigid  strip  footings  analyzed  In  their  study  showed  the  gradual 
development  of  the  plastic  zone,  the  corresponding  displacement  fields,  and 
tho  accumulation  of  footing  settlements  as  the  full  shearing  resistance  of  the 
foundation  material  was  mobilized.  The  load  settlement  rur.  js  approached 
asymptotically  the  bear'ng  capacity,  which  was  In  ayr  ’moot  with  ’he  ultimate 
load  predicted  by  plastic  theory. 


*  In  a  boundary  value  problem,  one  seek)  to  determine  the  response  (that  Is,  stresses, 
strains,  displacements)  of  a  system  (that  Is,  a  deformable  mass  of  soli)  to  a  specified  set 
of  boundary  conditions  (that  Is,  rpplled  loads  or  displacement). 


While  the  medium  properties  in  the  study  performed  by  Harper  and 
Ang  corresponded  to  steel  rather  than  soil,  and  therefore  the  conclusions  are 
based  on  the  behavior  of  such  an  idealized  "foundation"  material,  the  results 
do  provide  valuable  insights  into  the  behavior  attendant  on  the  loading  of 
soils,  The  lumped  parameter  modal  can  tolerate  any  type  of  stress-strain 
relationship,  and  future  studies  can  introcuce  stress-strain  proporties  snd 
yield  phenomena  better  representing  soli  materials, 

Although  the  basic  calculations  used  In  the  Harper  and  Ang  numerical 
method  are  very  simple,  the  large  number  of  calculations  required  to  solve 
even  a  simple  practice  problem  preclude  the  use  of  hand  calculations.  There¬ 
fore  the  entire  procedure  for  handling  plane  strain  problems  of  contained 
plastic  flow  in  an  elastic— perfectly  plastic  continuum  has  been  coded  in 
FORTRAN  II  for  use  on  the  IBM  7094  digital  computer.  The  computer 
program  used  In  this  study  Is  an  extensive  revision  of  a  program  written  for 
Whitman  (1964)  at  the  Stanford  Research  Institute.  By  revising  the  program, 
compilation  and  execution  time  hes  been  reduced,  some  logical  errors  have 
been  corrected,  and  the  output  has  been  reorganized  to  facilitate  its  use, 

This  appendix  presents  the  numerical  method  developed  by  Harper 
end  Ang  and  documents  tha  revised  Whitman  program,  No  attempt  has  been 
made  to  extend  the  theory, 


LUMPED  PARAMETER  MODEL 

Harper  and  Ang  (1963)  note  that  there  are  vary  few  solutions  for 
problems  of  contained  plastic  flow,  That  Is,  there  are  very  few  solutions 
which  trace  the  development  of  the  stress  and  strain  patterns  In  e  body  of 
material  from  the  time  that  yielding  first  develops  at  soma  point  In  the 
material  until  the  deformations  finally  Increase  beyond  all  bounds, 

One  criterion  that  has  been  used  In  tha  selection  of  a  mathematical 
model  Is  that  there  be  mathematical  consistency  between  the  finite  difference 
equations  governing  the  behavior  of  the  model  end  the  differential  equations 
governing  the  behavior  of  the  continuum.  This  means  that  the  equations  for 
stresses,  strains,  equilibrium,  and  compatibility,  which  are  derived  directly 
from  the  model,  should  be  the  seme  as  a  set  of  finite  difference  equations  of 
the  corresponding  differential  relations  governing  the  continuum,  Harper 
and  Ang  point  out  that  If  this  requirement  Is  met,  the  requirement  of  equal 
deformations  In  the  model  and  corresponding  continuum  will  be  automatically 
satisfied,  The  model  proposed  by  Harper  and  Ang,  and  used  herein,  satisfies 
this  criterion. 


49 


Description  of  the  Model 


Figure  C-1  shows  the  lumped  parameter  model  used  in  the  procedure 
developed  by  Fiarper  and  Ang,  This  model  consists  o'i  mass  points  and  stress 
points, 

The  mass  of  the  material  Is  assumed  to  be  concentrated  at  the  mass 
points  and  the  strains  within  the  material  arc  given  In  terms  of  the  displace¬ 
ments  of  the  mass  points,  Each  of  the  mass  points  is  connected  through  the 
stress  points  to  the  neighboring  mass  points,  Three  components  of  stress  and 
strain  are  defined  at  each  stress  point  (two  normal  components  and  a  shear 
component)  Displacements  In  the  continuum  are  defined  only  at  the  mass 
points  while  stresses  and  strains  are  defined  onl-  at  the  stress  points.  In  the 
model  representation,  springs  are  shown  at  the  stress  points  for  conceptual 
purposes  only,  The  stresses  at  each  stress  point  are  related  to  the  movements 
of  the  four  surrounding  mass  points  using  a  finite  difference  version  of  the 
stress-strain  laws  for  the  continuum,  No  attempt  is  made  to  define  a  set  of 
springs  whose  action  would  be  equivalent  to  the  action  of  tha  continuum. 

Actual  computations  are  carried  out  In  terms  of  the  forces  F„,  Fy, 

F,  at  each  stress  point  (Figure  C-2).  These  forces  are  equal  to  the  corre¬ 
sponding  components  of  the  stress  tensor  multiplied  by  the  appropriate 


Figure  C-1,  Lumped  perimeter  model  for  continuum. 


BO 


f 


upon  unrounding  m*w  point* 


eeord'.nitM 


¥  * 


Flgur*  C-2,  Foret*  at  itrtn  point  of  lumptd  ptrairtttr  modtl. 


Harper  and  Ang  and  Whitman  note  that  there  are  two  Important 
advantages  of  tho  modal  configuration,  First,  all  elements  of  the  strain  and 
stress  tensors  are  defined  at  tho  same  point,  This  Is  Bspeclally  Important  In 
extending  the  use  of  the  model  to  problems  of  plasticity,  Second,  the 
horizontal  and  vertical  boundaries  of  the  model  contain  only  mass  points. 
Thus,  boundary  conditions  given  In  terms  of  either  external  load"  acting  on 
the  mass  points  or  specified  displacements  of  the  mass  points  can  be  applied 
with  equal  ease, 

Relation  of  the  Model  to  the  Finite  Difference 
Equations  of  the  Continuum 

The  notation  and  example  presented  by  Harper  and  Ang  are  used  to 
Illustrate  the  relation  of  the  numerical  model  to  the  finite  difference  equations 
of  the  continuum.  The  following  notation  is  used; 

1,  Superscript  letters  refer  to  stress  point  locations, 

2,  Subscript  letters  x  and  y  refer  to  the  direction  of  the  exes. 

3,  Subscript  numbers  refer  to  mass  point  locations, 

4,  Displacement  components  In  the  x  and  y  directions  are  given  by 
u  and  v,  respectively. 

6.  Sign  convention  Is  that  shown  in  Figure  C-2. 
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For  Figure  C*1  the  components  of  the  strains  at  a  typical  stress  point,  a,  are 
defined  as  follows: 

UB4  -  U43 
-  - ! - 
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where  e  ■  normal  strain 
7  ■  shear  struln 

3  ■  diagonal  distance  between  mass  points 

These  strains,  which  are  derived  directly  from  tht  model,  are  Identical  to  the 
finite  difference  expressions  for  the  differential  uraln-dlsplacement  relations 
from  tha  classical  theory  for  plane  continue  under  small  deformations: 
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Using  the  sign  convention  of  Figure  C-2,  the  equation  of  equilibrium 
In  the  x  direction  for  a  typical  Interior  mass  point,  43,  is 


IF*  -  F«>  +  IF*  •  F*  )  +  -  0  (C-2) 

where  X  ■  body  force  per  unit  volume 

F  -  component  force  at  stress  point. 

The  volume  of  a  parallelepiped  of  unit  thickness  and  area  X3  ■  63/2  Is  corn 
sldered  concentrated  at  each  mass  point,  If  the  thickness  of  the  model  Is 
taken  as  unity  In  the  1  direction,  forces  at  the  stress  point  a  can  be  obtained 
from  the  stresses: 
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Using  Equations  C-2  and  C-3,  the  following  aquation  of  equilibrium,  In  terms 
of  stresses,  Is  obtained: 
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A  similar  equation  Is  obtained  for  the  y  direction: 
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These  equilibrium  equations,  C-4a  and  C-4b,  are  Identical  to  the 
finite  difference  expressions  for  the  differential  equations  of  equilibrium 
governing  the  corresponding  continuum: 
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Harper  and  Ang  state  that  "the  strains  In  the  model  will  necessarily 
satisfy  the  compatibility  relation,  since  strain  compatlDlllty  Is  essentially  a 
requirement  placed  on  the  three  components  of  strain  In  order  to  Insure  that 
they  correspond  to  a  physically  possible  displacement  configuration,  The 
model  deals  directly  with  displacements,  and  the  strains  are  defined  directly 
In  terms  of  these  displacements.  Hence,  It  can  be  expected  that  the  strains 
derived  from  the  displacements  of  the  model  will  exactly  satisfy  the  com¬ 
patibility  condition." 
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It  Is  also  possible  to  express  the  equations  of  equilibrium  in  terms  of 
displacements.  For  an  elastic  continuum,  the  stress-strain  relations  for  plane 
strain,  as  given  by  Prager  and  Hodge  (1951),  are 
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where  E  -  modulus  of  elasticity  of  the  material 
v  ■  Poisson's  ratio  of  the  matarlal 

Combining  Equations  C-1,  C-3,  and  C-5,  the  three  force  components  at  stress 
point  a  In  terms  of  displacements  are 
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The  three  expressions  In  Equation  C-6  are  Hooke's  stress-strain  relationships 
for  plane  strain  In  terms  of  displacements,  Substitution  of  these  rind  similar 
relations  for  the  forces  originating  at  the  other  stress  points  (b,  c,  and  d) 
into  Equation  C-2  results  In  the  following  equation  of  equilibrium  in  thex 
direction,  In  terms  of  displacement  as  given  by  Harper  and  Ang: 
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A  similar  equation  exists  for  equilibrium  In  the  y  direction,  Note  that 
Equation  C-7  is  the  same  as  a  finite  difference  equation  for  the  differentia1 
equation  of  equilibrium  governing  the  continuum: 


In  general,  boundary  conditions  (for  either  continue  or  discrete 
models)  can  be  of  two  types:  either  the  forces  acting  along  some  boundary 
<y  the  displacements  on  the  boundary  are  prescribed.  It  has  been  previously 
nored  that  the  model  Is  suited  to  either  type  of  prescribed  condition,  Some 
of  the  more  common  boundary  conditions  which  may  be  Imposed  on  the 
model  are  described  In  the  following  paragraphs, 

Symmatrlo.  If  the  continuum  Is  known  to  possess  symmetry  about  a 
vertical  axis  through  a  column  of  mass  points,  B3  shown  In  Figure  C-1‘,  then 
the  boundary  condition  on  the  right  edge  of  the  model  may  be  specified  as 
follows: 


uie  "  vu 

I  *  1,2, I*, ,6 

vie  ■  U|4 

Infinite.  The  model  may  be  used  to  simulate  an  Infinite  half-space. 

In  Figure  C-1,  consider  the  problem  of  Imposing  boundary  conditions  on  the 
fur-left  column  of  mass  points.  For  vertical  loadings  which  are  symmetric 
about  the  centerline,  It  has  been  ossumed  that  the  horizontal  displacements 
of  this  far-left  column  are  zero,  and  that  the  vertical  displacements  of  this 
far-left  column  will  be  equal  to  the  vertical  displacements  of  the  column  of 
mass  points  Immediately  to  the  right  of  this  boundary  column,  When  these 
vertical  and  horizontal  motions  are  resolved  Into  displacements  In  the  x  and  y 
directions,  the  boundary  conditions  become 

uii  ■  y  +  vu) 
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Specified  by  Loading  Arrangement.  If  an  external  load  Is  to  bo 
■xoplled  to  the  top  surface  of  the  continuum,  tlm  model  will  have  the 
appropriate  concentrated  loads  applied  to  the  top  row  of  mas9  points.  These 
concentrated  loads  may  approximate  a  distributed  load  or  represent  actual 
concentrated  loads. 

Fixed,  If  It  Is  desired  to  hold  the  base  of  the  continuum  fixed 
against  displacement,  the  displacement  components  of  the  bottom  row  of 
mass  points  am  simply  set  oqual  to  zero. 

These  examples  Indicate  tho  manner  In  which  boundary  conditions 
are  prescribed  for  the  model,  A  variety  of  practical,  significant  conditions 
can  be  conceived,  but  an  extensive  treatment  of  possible  boundary  condi¬ 
tions  Is  beyond  the  scope  of  this  study, 


CONSTITUTIVE  EQUATIONS  FROM  THE  THEORY  OF  PERFECTLY 
PLASTIC  HO  LI  OS 

Pragar  and  Hodge  (1061 )  have  provided  a  careful  definition  for  the 
perfectly  plastic  solid. 

Any  constitutive  relationship  of  the  theory  of  plasticity  may  be 
divided  Into  the  following  three  parts: 

1,  Stress-strain  relations  for  the  elastic  region 

2,  Yield  criterion  to  define  the  Initiation  of  yielding 

3,  Stress-strain  relations  for  the  ploBtlc  region 

These  three  major  divisions  of  tho  theory  will  ba  discussed  after  the  associated 
assumptions  and  limitations  era  listed,  and  after  a  net  of  notation  that  will  be 
useful  In  the  dlsaussion  of  the  theory  Is  Introduced, 

Harper  and  Ang  (1063)  state:  "There  are  three  main  assumptions 
underlying  the  theory  of  perfectly  plasti :  material  used  In  this  Investigation, 
These  can  be  stated  ns  follows: 

"  1 ,  1 1  Is  assumed  that  tho  Mlses-Hencky  yield  condition  accurately 
determines  the  beginning  of  yield,  General  considerations  of 
Isotropy  end  symmetry  can  furnish  only  the  general  form  of  the 
yield  condition,  Beyond  this,  any  yield  condition  Is  a  hypothesis 
which  only  tests  can  justify. 

"2,  It  Is  assumed  that  there  Is  no  permanent  volume  change,  This 
assumption,  justified  on  the  basis  of  experimental  evidence  for 
metals,  leads  to  the  result  that  the  plastic  strain  Is  equal  to  the 
plastic  deviator  strain, 


"3.  During  the  plastic  flow,  It  l» 
assumed  that  the  devhtor 
strain  rate  tenujr  Is  proportion' 
al  to  the  Instantaneous  deviator 
stress,  This  Is  the  familiar 
Prandtl-Reuss  postulate, 

*  "In  addition  to  those  threa  main 
assumptions,  It  Is  possible  to  list 
several  other  restrictions  on  the 
theory: 

"4,  The  material  must  be  Isotropic, 
This  condition  Is  uaad  In  dovel* 
figure  C*9.  Stre»itraln  eurve  for  oping  the  general  form  of  the 

elaitle-psrfeotlv  pintle  materiel  vuw  eonc{|tlon, 

In  simple  tension  or  compression, 

"8,  There  le  no  work  hardening, 
and  the  material  followi  tha 
etrosutreln  diagram  of  [Figure  C*3]  when  subjaotad  to  simple 
tension  or  comnresilon, 

"0,  No  unloading  occurs,  Once  a  stress  point  has  yielded,  It  remains 
yielded  under  successive  Increments  of  external  losd, 

"7,  Tima  effects  of  loading,  such  as  creep,  are  Ignored, 

"8,  Displacements  ere  smell  so  that  the  smell  dsformitlon  theory  of 
elasticity  applies, " 

The  following  definitions  and  notation  are  Introduced  for  tha  purposa 
of  describing  the  pertinent  constitutive  equations  used  In  this  study: 
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Spherical  stress  tensor 
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Doviator  stress  tensor  =  SD 


«x 

Txy 

I 

*xy 

*v 

TV* 

Ty« 

*1 

where  »  ■r  mean  normal  stress  ■  1/3  (aK  +  <7V  +  ot) 

»x  ■  normal  x  component  of  S°  ■  e„  -  s 

*y  *  normal  y  component  of  S°  -  cv  -  s 

8,  *  normal  z  component  of  SD  <■  c,  -  s 


With  this  notation 


ST  *  S8  +  SD 


•„  +  »y  +  3,  ■  cx  +  oy  +  o,  -  3 1  ■  0 

Principal  normal  stresses  are  designated  by  cr,,  o2,  and  o3,  Principal  normal 
components  of  the  stress  deviator  sre 
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A  completely  similar  notation  exists  for  strains: 


Total  strain  tensor  “  ET 
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Devisor  strain  tensor 
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where 


mean  normal  strain  “  ‘  1/3  +  ey  +  e,) 

normal  x  component  of  E°  ■  eK  -  e 

normal  y  component  of  E°  ■  ev  -  e 

normal  z  component  of  E°  »  e,  -  e 


With  this  notation 


ET  ■  Es  +  E° 

e*  +  ev  +  e,  «  e„  +  ey  +  e,  -  3  a  ■  0 

Principal  normal  strains  are  designated  as  e, ,  e3,  and  e3.  Principal  normal 
components  of  the  strain  deviator  are 

e,  «  e,  -  e 

ea  ■  (j  -  j 

e3  -  a3  -  e 


Elastic  Stress-Strain  Relations 

In  the  elastic  range  the  relationship  between  the  elements  of  the  stress 
and  strain  tensors  Is  assumed  to  be  that  of  Hooke's  law.  It  Is  convenient  to 
express  this  linear  relationship  In  terms  of  the  elements  of  the  deviator  stress 
and  deviator  strain  tensors: 
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where  G 
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The  relationship*;  in  Equation  C-8  can  be  expressed  more  rnnrisply  as 

SD  -  2  G  tD  (C- 10) 

Note  that  the  expressions  In  Equation  C-8  or  Equation  C-10  are  not  six 
independent  relations  since  addition  of  ex  +  ey  +  »,  -  0  gives  an  identity. 
Therefore,  Equation  C-9  is  needed  to  give  a  complete  statement  of  Hooke's 
law. 

Yield  Criterion 

The  yield  point  is  determined  by  the  Mlses-Hencky  yield  criterion: 

Ja  ■  ka  (C-1D 

where  Jj  ■  second  Invariant  of  the  stress  deviator  tensor 
k  ■  yield  stress  In  simple  shear 

J3  is  defined  by 
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For  plane  strain  problems,  the  yield  condition  reduces  to 


(C-12) 
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Combining  fcquat.ons  C-1 2  and  C-y,  we  obtain  for  any  yielded  stress  point 


4 


Still  another  useful  form  for  the  yield  criterion  is  obtained  by  substituting  the 
expressions  in  Equation  C-5  into  Equation  C-1 2: 


<•*  ev>2  +  V  “ 

Using  Equations  C-1  for  the  strains  at  a  typical  stress  point,  a,  the  yield 
criterion  becomes: 


(UB4  -  U43  -  V53  +  v44}2  +  (uB3  -  tl44  +  VB4  -  V43)2 


(C-1 3) 


Plastic  Stress-Strain  Relations 

In  order  to  relate  stress  and  strain  In  a  material  which  is  undergoing 
plastic  flow,  It  Is  convenient  to  express  the  strain  tensor  In  terms  of  elastic 
and  plastic  components.  Single  primes  will  be  used  to  denote  an  elastic 
component,  and  double  primes  will  denote  a  p'astic  component,  Dots  will 
denote  rate  of  change  with  respect  to  increment  of  external  load, 

The  essential  nature  of  the  relations  between  stress  and  strain  during 
plastic  flow  is  given  by  Equations  C-14  and  C- 1 5.  The  assumption  of  no 
permanent  change  of  volume  Is  stated  mathematically  as 

e"  -  ~  (e"  +  e"  +  e'')  -  0  (C-1 4) 

This  Implies  that  the  plastic  strain  deviation  is  identical  to  the  plastic  strain, 
or 


The  assumption  that  during  plastic  flow  the  deviator  strain  rate  tensor 
Is  proportional  to  the  instantaneous  deviator  stress  tensor  is  stated  mathemat¬ 
ical  iy  as 
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(C-15) 
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where  $  ■  a  proportionality  factor 

The  expressions  in  Equation  C-15  are  in  the  same  form  as  the  elastic 
stress-strain  relations  given  in  Equation  C-8. 

The  basic  relationships  which  are  assumed  during  plastic  flow  have 
been  presented,  At  this  point  it  is  necessary  to  apply  these  relations,  along 
with  the  yield  criterion  (Equation  C-1 1 )  and  the  elastic  relations  (Equations 
C-8  and  C-9),  in  order  to  develop  the  final  relationships  between  the  stress 
rates  (incremental  stresses),  strain  rates  (incremental  strains),  and  instantaneous 
stresses. 

The  plastic  strain  rates  have  been  expressed  In  terms  of  stresses  by 
Equation  C-15.  Similarly,  the  elastic  strain  rates  are  expressed  in  terms  of 
stress  rates  by  differentiating  the  expressions  in  Equation  C-8  with  respect  to 
external  load: 
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Combining  the  elastic  and  plastic  strain  rates  gives  the  total  strain  rate: 
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(C-1 6) 


Note  that  these  relations  apply  only  during  plastic  flow,  that  Is,  when 
J2  “  k2  and  Ja  ■  0 
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In  order  to  eliminate  the  proportionality  factor  0  from  the  expressions 
in  Equation  C-16,  it  is  convenient  to  introduce  the  notation 

W  -  sx ex  +  svey  +  s,e,  +  Txy7xy  +  rx,7x,  h-  ry,7vl 

where  W  may  be  interpreted  as  the  rate  at  which  stresses  do  work  during  a 
change  of  shape,  and  to  note  that 

^2  -  M*  +  sysy  +8,3,  +  2rxyTxy  +  2 rx, txi  +  2tv,tv, 

By  multiplying  the  first  three  expressions  of  Equation  C-16  by  sx ,  sy ,  8,,  and 
the  last  three  by  2  rxy,  2  rxl,  2  ryI,  respectively,  and  adding,  we  get 

2GW  ■  «„•„  +  0ex2  +  «yiy  +  08y2  +  8,8,  +  08, 2 

+  2TxyTXy  +  2  0T„y2  +  2  rx,  r„,  +  2  0txi  , 

+  2ryIrvl  20fy,2 

-  jj  +  0(8„2  +  8y2  +  8, 2  +  2  fay2  +  2r„,2  +  2ry,a) 

*  +  2  (p 

But  during  plastic  flow,  Ja  ■  k2  and  j2  ■  0,  Hence,  2GW  ■  2 0k2  and 
0  -  G  W/k2. 

Substituting  this  value  of  0  Into  the  expressions  In  Equation  C-16,  It 
Is  possible  to  solve  for  the  deviator  stress  rates,  which  gives 

'  2G('"  '2k3’”)  '-V  ■  °  (f«  -  7  r») 

<  '  2G(*»  "  Q(1'"  '7’")  lc',7) 

'■  '  2°(,‘  -  *■)  '  °(fv  "p  rv) 
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To  obtain  the  total  stress  rates  it  is  necessary  to  add  the  deviator  stress 
rales  from  the  expressions  in  Equation  C-17  to  the  spherical  stress  rate,  which 
can  be  obtained  by  differentiating  Equation  C-9  with  respect  to  external  load: 


s  “  3Ke 

Adding  Equations  C-17  and  C- 1 6  results  in  the  total  stress  rates: 
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(C-10) 


The  expressions  in  Equation  C-19  give  the  desired  relationships  between  the 
stress  rates,  strain  rates,  and  Instantaneous  strosses. 

To  apply  the  expressions  in  Equation  C-19  to  the  numerical  model,  It 
Is  necessary  to  reduce  them  to  an  Incremental  form,  Note  that  for  plane  strain 
problems  the  number  of  relations  Is  reduced  from  six  to  three,  Therefore 


Aax  ■  Asx  +  As 

Affy  ■  Aly  +  As 

Ar  .  q  Ur  r\ 

UT*i  i  ar¥y  ^  T«yl 


(C-20) 
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For  plane  strain  onditions,  the  expressions  in  Equation  C-17  are  reduced  to 
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Ar„v  -  G  { A7 


/a  AW  \ 
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and  Equation  C-18  becomes 

As  ■  3KAe  ■  K(Ae„  +  Aev) 

The  increment  W  becomes 

AW  ■  SaAe*  +  sy  At*y  +  i,  Ae,  +  T^Ay^  (C-21) 
But  s,  ■  o,  -  y  +  oy  +  0,) 

Where,  for  plane  strain 

0,  -  K(0„  +  Oy) 

and  during  plastic  flow 
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Similar  expressions  are  obtained  for  Aoy  and  A txv: 


(C‘26) 

Combining  Equations  C-1  and  C-3  with  Equations  C-24,  C-26,  and 
C-26,  and  rearranging,  tha  thraa  force  components  at  the  stress  point  a  In 
terms  of  olsplacennanta  are 

AFk  ■  A(Auj4  -  Au43)  +  B(Av33  -  Av44) 

-  C(Aum  -  Au44  +  Av#4  -  Av43)  (C-27) 

AFy  ■  B(AuI4  -  Au43)  +  A(AV|3  -  Av44) 

+  C(Au4j  -  Au44  +  Av34  -  Av43)  (C>28) 

APliV  -  -  C(AuB4  -  Au43)  +  C(Av#j  -  Av44| 

+  D(AU|,  -  Au44  +  Av,4  -  Av43)  (C-20) 
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v;here  the  A  denotes  change  during  an  external  load  increment  and 
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Equations  C-27,  C-20 .  and  C-29  era  the  relationships  with  which  the 
Incremental  force  components  In  a  plast-c  region  are  computed,  These 
Incremental  force  components  are  added  to  the  existing  force  components 
(F„,  Fv,  and  F„v)  at  a  yielded  stress  point  to  obtain  the  total  forces. 

In  order  to  compute  the  quantities  Au  and  Av  which  appear  In 
Equations  C-27,  C  2B,  and  C-29,  two  sets  of  displacements  corresponding  to 
two  consecutive  load  levels  are  required.  One  set  of  displacements  Is  the  set 
which  Is  being  generatod  for  the  current  load;  the  other  set  Is  that  computed 
for  the  previous  load.  The  quantities  Au  and  Av  are  computed  as  the 
difference  In  displacements  determined  for  these  two  loads, 

Equations  C-27,  C-28,  and  C-29  are  linear:  that  Is,  If  A,  B,  C,  and  D 
are  constants,  then  the  changes  In  force  are  linear  with  regard  to  changes  In 
displacement,  Actually,  the  equations  ere  not  linear  since  In  general  F„,  Pyi 
FKy  will  change  during  plastic  flow.  Thus,  It  Is  necessary  to  proceed  In  a 
step-wise  linear  fashion. 


General  Form  of  the  Pleatlulty  Equations 

In  general,  the  application  of  the  plasticity  equations  to  the  Harper 
and  Ang  model  Is  closely  associated  with  the  three  stages  of  materiel  behavior 
presented  In  the  previous  sections.  The  plasticity  relations  which  describe 
those  three  Hugos  of  material  behavior  In  terms  of  stress  point  forces  and 
mass  point  dlsplacamants  are  used  \o  describe  the  behavior  of  the  numerical 
modal. 
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Thusn  relations.  which  have  been  developed  in  the  previous  sections, 
can  be  generalized  by  adopting  the  folln  wing  subscript  notation: 

l.R  ■  lower  right 
UL  ■  upper  left 
LL  ■  lowor  left 
UR  ■  upper  right 

Forca- Displacement  Relation*:  Elastic  Caie.  The  expressions 
(Equations  C-6)  which  relate  olastic  stress  point  forces  to  mass  point  dis¬ 
placements  can  be  generalized  as 

F»  "  arr+7)“aTi  [(1  “ u){u^  *  1  ,>(v^  ‘  vu«>] 

FV  '■  2(rv-kr-  a'y)  [(1  “  M,(v“  ■  vu*)  +  ‘'K*  ■  Uuu>]  <C  30> 

P"v  "  4(1  A  mS  ^u.  '  uu*  vtR  fc  vuu> 

Yield  Criterion.  The  yield  criterion  (Equation  C- 13)  In  termiof  men 
point  displacement!  can  be  generalised  m 

(uUB  -  uUL  -  VuL  +  Vg*)* 

+  (Utk  “  Uu«  +  vLB  -  vuu)«  ■  (C'31) 

Force-Displacement  Relitlonn  Pleitlo  Ceie.  The  expressions 

(Equations  C*I47,  C-2B,  end  C-20)  which  relate  Inorementil  pleitlo  stress  point 
forces  to  mess  point  displacements  csn  be  generallisd  as 

AF(  *  A(AUgR  -  AUgg)  +  IfAVn,  ■  AVg r  ) 

“  C(AUgg  ■  AUyp  ■  AVy^  “■  AVgg) 

AFy  ■  B(AuL„  -  Auuv,}  +  A(Avu  -  Avy„) 

+  C(Augy  ■■  Auy^  »  Avgn  -  Avgg)  1C  32) 

AF,y  •*  •  CUUiR  «  Augg)  ♦  CfAVgj,  •  AVg^) 

+  D(AUU  -  Augq  •  AVgn  «  AVyi.) 
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Equations  C-30,  C*3 1 ,  and  C-32  are  the  form  of  the  plasticity  relations 
actually  used  In  the  numerical  method, 


NUMERICAL  METHOD 

When  a  problem  In  continuum  mechanics  Is  replaced  by  a  corresponding 
problem  In  particle  mechanics  Involving  a  discrete  mofel,  the  displacements  to 
satisfy  equilibrium  must  be  determined,  One  method  is  to  write  and  solve  the 
set  of  simultaneous  linear  algebraic  equations  (equations  similar  to  Equation 
C-7)  for  the  unknown  displacement  components  u  and  \  at  each  mass  point, 
HBrper  and  Ang  (1963)  note  that  such  an  approach  has  significant  disadvan¬ 
tages,  The  preparation  of  the  equations,  whether  It  Is  done  by  hand  or  by  on 
Intricate  program  for  the  computer,  Involves  a  considerable  amount  of  labor, 

In  addition,  oven  with  machines  as  large  as  the  IBM  7094,  the  number  of 
actuations  which  can  be  solved  by  the  standard  library  subroutines  Is  limited, 
Perhaps  nost  Important,  the  changes  In  t m  coefficients  for  the  displacements 
resulting  from  the  yielding  of  one  or  rnora  stress  points  are  not  easy  to  deter¬ 
mine. 

A  more  flexible  and  practical  approach  to  the  problem  as  suggustod  by 
Harper  and  Ang  (1963)  and  Whitman  (1964)  Is  to  employ  a  relaxation 
technique,  Such  an  approach  eliminates  completely  the  preparation  of 
simultaneous  equations,  and  can  handle  several  thousand  displacement  com¬ 
ponents.  An  additional  advantage  of  the  relaxation  method  Is  the  physical 
meaning  that  can  be  attached  to  each  stop  of  the  procedure.  This  Is  very 
helpful  In  determining  plastic  forces  and  displacements. 
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Relaxation  Procedure 


The  relaxation  procedure  used  for  determining  the  equilibrium 
displacements  is  summarized  by  the  flow  diagram  in  Figure  C-4.  Following 
the  discussion  given  by  Harper  and  Ang,  all  mass  points  of  the  model  are 
initially  in  equilibrium  with  zero  displacements  and  no  external  load.  The 
first  increment  of  external  load  is  then  applied  to  any  or  all  moss  points, 
thus  destroying  the  equilibrium  of  the  loaded  mass  points.  The  following 
operations  are  then  performed  for  each  mass  point  of  the  model. 

The  forces  acting  on  a  mass  point  are  determined  as  follows: 

External  forces  acting  on  the  mass  point  are  given  as  a  port  of  the  loading 
pattern  applied  to  the  model.  Internal  forces,  originating  at  the  stress  points, 
are  determined  uniquely  by  Equation  C-30  in  the  elastic  range  from  the 
displacements  surrounding  the  stress  points,  After  a  stress  point  has  yielded, 
the  force  components  at  that  stress  point  are  determined  both  by  the 
surrounding  displacements  and  the  past  history  of  that  particular  stress  point. 
Incremental  plastic  forces  (Equation  C-32)  are  then  added  to  the  last  set  of 
equilibrium  forces  at  the  stress  point  to  obtain  the  current  total  plastic  forces 
acting  at  the  yielded  stress  point, 

After  the  forces  acting  on  a  given  mass  point  are  determined,  a 
summation  of  all  the  forces  acting  In  thex  direction  Is  mode,  In  general, 
this  will  result  In  a  residual  force  which  is  an  Indication  of  the  amount  by 
which  the  mass  point  Is  out  of  equilibrium  In  tho  x  direction.  The  mass  point 
Is  then  displaced  through  a  small  dlstanco  In  the  x  direction  equal  to  the 
product  of  the  residual  force  and  a  flexloillty  coefficient, 

Similar  operations  are  performed  for  tho  y  direction,  These  operations 
place  tho  current  mass  point  in  equilibrium,  though  in  general  the  equilibrium 
of  surrounding  mass  point.,  will  bo  upset  by  a  small  amount,  The  procedure  Is 
repaatod  for  each  moss  point  until  every  mass  point  has  been  moved  once  in 
the  x  direction  and  once  In  the  y  direction,  thus  completing  one  cycle  of 
relaxation, 

After  every  relaxation  cycle,  each  muss  point  la  Inspected  to  determine 
if  It  Is  in  equilibrium.  If  not,  the  relaxation  process  is  repeated  until  all  mass 
points  are  In  equilibrium  within  the  accuracy  prescribed  by  a  convergence 
criterion.  After  all  moss  points  arc  In  equilibrium,  all  the  stress  points  am 
inspected  for  yielding  by  Wises- Hencky  yield  criterion  (Equation  C-3 1 )  and 
tho  yielded  regions  are  recorded,  All  the  displacements  anti  forces  for  the 
equilibrium  configuration  |ust  obtained  are  also  rocorded. 

If  desired,  the  external  load  is  given  a  now  increment,  and  the 
complete  procedure  Is  repeated  for  each  load  increment  In  order  to  trace 
the  development  of  plastic  yielding  from  one  stress  point  to  another, 


71 


1 

i 


I 


1 


Flour#  C-4.  Flow  diagram  for  relaxation  procedure. 
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This  relaxation  procedure  has  been  coded  for  use  on  the  digital 
computer.  Although  the  basic  calculations  are  sirnpie,  a  complex  program  is 
needed  to  handle  all  of  the  logical  decisions  that  the  computer  must  make  as 
it  distinguishes  between  yielded  and  unyielded  stress  points  and  computes 
the  magnitude  of  the  load  increment  needed  to  just  cause  yielding  at  one 
additional  stress  point. 

Program  PEHFPLAS  II 

The  computer  program  used  in  this  study  is  a  revision  of  the  program 
(PERFPLAS  I)  developed  for  Whitman  (1964)  at  the  Stanford  Research 
Institute.  The  revised  program  has  been  named  PERFPLAS  II  (PER Factly 
PLAStic)  and  Is  written  In  FORTRAN  II  for  the  IBM  7094  digital  computer. 

At  present,  program  PERFPLAS  II  will  permit  the  use  of  up  to  a 
36  x  36  grid  of  mass  points  (36  x  35  stress  points).  The  allowable  grid  size 
can  be  Increased  by  Increasing  the  storage  allocated  by  the  program's 
dimension  statements, 

The  boundary  conditions  can  be  of  four  types:  free,  fixed,  reflected, 
or  Infinite.  A  free  boundary  Is  one  in  which  the  mess  points  can  move  In  any 
direction  and  a  fixed  boundary  is  one  In  which  the  mass  points  cannot  move. 

A  reflected  boundary  Is  a  line  of  symmetry  between  two  halves  of  a  symmet¬ 
rical  problem  for  which  It  is  necessary  to  solve  only  one  hulf,  An  Infinite 
boundary  Is  an  approximation  of  the  condition  far  from  the  loaded  area,  In 
that  mas3  points  can  move  only  parallel  to  the  boundary  in  such  a  manner 
that  normal  lines  remain  normal  to  the  boundary, 

In  PERFPLAS  II,  the  left  boundary  is  limited  to  a  fixed  or  an  Infinite 
boundary.  The  right  boundary  Is  limited  to  a  fixed,  a  symmetrical,  or  an 
Infinite  boundary.  The  top  boundary  Is  United  to  a  free  boundary, 

The  basic  step  In  the  program  Is  the  application  of  a  trial  load  Increment 
and  the  successive  relaxation  of  mass  points  until  the  system  comes  to  equili¬ 
brium.  The  program  keeps  track  of  stress  points  which  have  yielded  and  uses 
the  plastic  equations  accordingly  (Equation  C-32) . 

Two  load  options  have  been  programmed  Into  PERFPLAS  II,  These 
load  options  are  the  same  when  the  material  deformations  are  elastic,  arid 
differ  only  after  the  material  has  yielded  at  some  point.  Initially,  a  small  load 
(known  to  be  less  than  the  load  required  to  cause  first  yielding)  Is  applied, 

When  the  specified  standard  load  Increment  is  used,  the  initial  load  Is  Incre¬ 
mented  until  the  first  stress  point  yields,  Thereafter,  the  two  load  options 
differ; 


1.  When  the  first  load  option  Is  used,  the  load  will  continue  to  bo 
Incremented  by  the  specified  standard  load  Increment  until  the 
load  attains  the  specified  maximum  value, 


73 


2.  W hnn  the  second  toad  option  is  used,  'he  load  win  incremented 
by  the  specified  standard  load  increment,  but  each  time  a  new 
stress  point  yields  the  load  will  be  reduced  to  a  level  which  will 
just  cause  this  stress  point  to  yield,  This  process  continues  until 
the  load  attains  or  slightly  exceeds  the  specified  maximum  load. 

The  relations  required  to  downgrade  the  load  when  the  second  load  option  i; 
is  used  are  presented  in  Appendix  D 

Whenever  an  adjustment  is  to  be  made  in  the  results  of  a  trial  load 
increment,  the  program  determines  the  amount  by  which  tho  load  increment 
must  be  changed  and  then  adjusts  all  computed  displacements  and  forces 
■accordingly,  By  using  the  second  load  option,  the  propagation  of  yielding 
within  the  material  may  be  studied  In  detail  (one  stress  point  at  a  time). 
However,  the  execution  time  of  the  program  will  be  Increased.  It  was  found 
convenient  to  use  the  first  load  option  for  the  first  run  of  all  new  problems, 

Iri  the  range  where  plastic  deformations  occur,  the  load  Increments 
should  be  sufficiently  small  so  that  the  development  of  the  plastic  zone  Is 
gradual.  This  keeps  the  errors  due  to  linear  extrapolation  to  a  minimum, 

As  mentioned  previously,  approximate  linearized  equations 
(Equation  C-32)  are  used  to  compute  the  forces  at  already  yielded  stress 
points.  For  this  reason,  the  yield  condition  may  be  exceeded  at  these  stress 
points  at  the  end  of  an  Increment  of  plastic  straining,  and  a  correction  must 
be  Introduced  to  bring  the  forces  back  to  the  yield  surface,  Tho  relations 
required  to  make  this  correction  are  presented  In  Appendix  D.  In 
RERFPLAS  II,  subroutine  CORRECT  makes  this  correction,  Subroutine 
CORRECT  Is  applied  after  equilibrium  has  been  established  under  the  given 
load  Increment.  This  subroutine  alters  the  forces  at  each  previously  yielded 
stress  point  without  altering  any  displacements,  As  a  result,  the  surrounding 
mass  points  are  slightly  out  of  equilibrium. 

A  convergence  criterion  Is  needed  to  specify  what  accuracy  Is  desired 
In  the  Iteration  process  following  each  applied  load  increment,  The  criterion 
adopted  Is  that  the  ratio  of  the  next  Incremental  adluotment  In  displacement 
to  the  total  displacement  at  the  mess  point  should  be  loss  than  or  equal  to  a 
prescribed  value,  The  next  incremental  adjustment  In  displacement  is 
determined  by  multiplying  the  unbalanced  force  acting  on  a  mass  point  by  an 
elastic  flexibility  coefficient,  Displacements  are  adjusted  In  both  the  x  and  y 
directions  to  bring  the  mass  points  into  equilibrium,  and  the  convergence 
criterion  Is  applied  to  all  mass  points, 

The  elastic  flexibility  coefficient  for  a  mass  point  Is  developed 
subsequently,  This  coefficient  Is  used  for  both  elastic  and  plastic  deforma¬ 
tions,  As  a  result,  the  number  of  Iterations  required  to  establish  equilibrium 
Increases  sharply  after  a  large  number  of  3tress  points  have  yielded. 
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In  program  PERFPLAS  II,  the 
mass  points  are  identified  by  a  double 
subscript  (I, J).  I  indicates  the  vertical 
position  of  the  mass  point  (the  row) 
and  increases  vertically  downward, 

J  indicates  the  horizontal  position  of 
the  mass  point  (the  column)  and 
increases  from  left  to  right,  Stress 
points  are  identified  by  a  double  sub¬ 
script  (l,J)  corresponding  to  the  mass 
point  above  and  to  the  left  of  the  stress 
point.  This  subscript  notation  Is 
Illustrated  In  Figure  C-6, 

Using  Program  PERFPLAS  I, 

Whitman  and  Hoeg  (1965)  have 
studied  the  performance  end  accuracy 
of  the  Harper  and  Ang  mathematical 
model.  From  this  study,  it  is  suggested  that  a  convergence  criterion  of  10'B 
or  10*8  should  be  used,  Christian  (1965'  has  shown  that  the  number  of 
Iterations  required  for  equilibrium  Increases  sharply  If  Poisson's  ratio  exceeds 
0,45.  This  Is  because  the  expression  (1  -  2v)  appears  In  the  denominator  of 
the  elastic  force  displacement  equations  (Equation  C-30).  As  Poisson's  ratio 
approaches  0,6,  the  expression  (1  -  2  v)  approaches  zero, 

The  FORTRAN  deck  of  program  PERFPLAS  II  can  be  compiled  In 
less  than  3  minutes  on  the  IBM  7094,  Compiling  time  nan  be  eliminated  by 
using  a  binary  deck  on  production  runs,  Execution  time  Increases  (for  e 
given  value  of  Poisson's  ratio)  as  the  number  of  mass  points  Increases  and  as 
the  number  of  yielded  stress  points  Increases,  A  detailed  study  of  execution 
time  has  not  been  made. 

The  user  of  program  PERFPLAS  M  must  provide  certain  data 
describing  the  problem,  The  data  can  be  In  any  consistent  set  of  units,  and 
the  output  will  be  In  the  same  units,  A  date  InDut  guide  is  given  In 
Appendix  E, 

A  glossary  of  notation  used  In  PERFPLAS  II  Is  given  In  Appendix  F 
and  a  listing  of  the  FORTRAN  deck  Is  given  In  Appendix  G,  Also  listed  In 
Appendix  G  Is  the  Input  data  for  the  sample  problem  which  will  be  discussed 
subsequently.  Program  output  will  be  discussed  In  connection  with  the 
sample  problem. 


75 


!  SAMPLE  PROBLEM 

A  sample  problem  has  been  selected  to  familiarly,-  the  reader  with 
the  numerical  method  and  program  PERFPL.AS  II. 

The  problom  selected  Is  the  seme  bs  that  presented  by  Harper  and 
■  Ang  (1903).  This  elementary  problem  (Figure  C  6)  was  selected  because  it 

f  can  be  readily  solved  by  hand.  Therefore,  the  numerical  procedure  can  be 

1  Illustrated  in  detail. 

j  In  Figure  C-6,  only  mass  points  12  and  22  are  free  to  move,  and  due 

1  to  symmetry  about  a  vortical  line  through  these  mass  points,  the  u  and  v 

displacements  at  a  mass  point  are  equal: 
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u12  =  V12  U22  ~  v22 

Hence  there  are  only  two  unknown  displacements,  u12  and  u22,  By  means  of 
the  material  constants,  dimensions,  and  loading  shown  in  Figure  C-6,  it  is 
possible  10  write  two  simultaneous  linoar  algebraic  equations  (similar  to 
Equation  C-7)  for  the  elastic  behavior  of  the  system  in  terms  of  the  two 
un'.nowns,  U!  2  and  u22.  Solution  of  these  two  equations  yields 

u12  a  v22  ■  1.429  x  1 0'2  inches 

(C-33) 

u22  ■  v22  -  3,571  x  10*3  inches 

These  values  will  now  be  used  to  measure  the  progress  of  the 
relaxation  procedure,  Converting  these  displacements  to  elastic  force  com¬ 
ponents  at  stress  point  12  by  using  Equation  C-30  gives 

(F«),a  -  -7.857  klpi 

(Fv)u  -  -0.714 kip  (C-34) 

(F„v)ia  -  -2.143  kips 

Before  the  systematic  relaxation  procedure  Is  begun,  It  Is  first 
necessary  to  convert  external  loads  to  concentrated  loads  for  application  at 
the  loaded  mass  points  and  to  determine  the  flexibility  coefficients  for  each 
mass  point,  For  example,  If  an  external  vertical  pressure  of  14.14  ksl  Is 
acting  on  the  top  surface  of  the  modal  shown  In  Figure  C-6,  the  concentrated 
vertical  force  acting  on  mass  point  1 2,  which  arises  from  tills  pressure  acting 
over  a  distance  of  X/2  *  1/2  Inch  on  either  side  of  mass  point  12,  is 

Pv  -  (14.14  ksl)^  +  j  Inch^  (1  Inch)  -  1 4.1 4  kips 

where  the  thickness  of  the  model  In  1  inch.  This  vertical  force  is  then 
resolved  into  components  in  the  x  and  y  directions  for  application  to  mass 
point  12: 
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The  flexibility  coefficient  for  a  mass  point  is  obtained  by  determining 
the  deflection  caused  at  the  mass  point  by  3  un>t  component  force.  All 
adjacent  mass  points  are  assumed  fixed,  and  the  flexibility  coefficient 
determined  will  be  for  the  direction  of  the  applied  unit  force,  For  exnmDlo, 
a  unit  force  of  1  kip  applied  in  the  x  direction  at  mass  point  22  is  resisted  by 
internal  force  components  at  stress  points  11,  12,21,  and  22,  By  using  the 
sign  convention  shown  in  Figure  C-2  and  by  summing  the  forces  In  the  x 
direction  at  muss  point  22,  we  get 


(IF, 


22 


-  (F. 


'22 


-  (F, 


1 1 


(F 


*V'2i 


-  (F 


xv'u 


+  1 


Expressing  (F„)  ,  (F„)  ,  <P»,V)21 ,  and  (F,v)ia  in  terms  of  displacements  by 
means  of  Equation  C-30,  and  noting  that  all  displacement  components  except 
u22  are  zero,  we  get 


or 


u22 


.  r  e(i  - v)  +  _jj _ _ 

[d  +  w<1  -  2u)  4(1  +  P)JU2: 

2(1  +  m)(1  -  Iv)  .  ,f  ,  ,, 

"  . 13""-  4mT1 .  "  (f"!22  "  <fV 


(C-35) 


Equation  C-36  gives  ihe  flexibility  coefficient  for  a  typical  intorlor  mass 
point  In  the  x  direction  and  has  units  of  inches/klp.  Due  to  symmetry,  this 
Is  tilso  the  flexibility  coefficient  for  a  typical  interior  moss  point  in  the  y 
direction, 

From  Figure  C-1  it  can  be  seen  that  the  flexibility  coefficient  for 
mass  point  12,  which  is  a  typical  top  boundary  mass  point,  is  twice  that  of 
the  flexibility  coefficient  for  mass  point  22,  or 


12 


2(\ 


22 


4(1  +  t>)(1  -  2 v) 
(3  -  4  k)E 


If  E  and  v  take  on  the  values  1,000  ksl  and  0,26,  respectively,  as 
shown  In  Figure  C-u,  these  flexibility  coefficients  become 


If*),,  -  <fy)ia  •  1.25  x  10-3  In./klp 


l\>„  -  (V„  -  6.25  x  10*4  In./klp 


With  the  values  for  the  concentrated  external  loads  and  the  flexibility 
coefficients  known,  It  Is  possible  to  begin  the  relaxation  procedure,  The 
Allowing  step  numbers  refer  to  the  flow  diagram  of  Flgurp  C-4, 
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Qporatlon 


Sotu12  *  t a  ■  u22  *  v22  ■  0,  AIsosbi  forco  components  -  0. 

Apply  Increment  of  nxtormi!  loud  to  muss  point  12 

-  10  kips 

(Py)1}  -  lOklpi 
Begin  with  muss  point  12. 

No  stress  point  has  yielded,  since  all  stress  components  ere 
Initially  zero.  Goto  6b, 

On  the  first  cycle,  all  forco  components  ore  loro  since  the  mass 
points  have  not  moved, 

On  tho  first  cycle,  only  external  forces  aro  nonzero,  Her.uti 

IP.),,  -  IP.),,  -  +10 klpv 

IPV)„  "  IPy)„  "  +10klp« 

Nowuu  -  old  u, a  +  (f*),a S;iFM)ia 

uia  -  0  +  (0.00120)00)  -  0,0128  Inch 
Similarly 

v„  *  0  +  (0.00125)110)  «  0,0128  Inch 

Note  that  those  displacements  of  muss  point  12  destroy  the 
tujiillltirliiin  of  mess  point  22. 

Tho  current  moss  point,  12,  Is  not  tho  hist  muxs  point,  Go  to  9, 
Toko  muss  point  22,  Go  to  4, 

Again  no  stress  point  luts  yielded,  since  yielding  can  occur  only 
alter  eqnlllbi  lum  Is  reuclvid,  Go  to  bl>, 
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Operation 


Foccg  components  at  stress  points  11  nnrl  12  are  compntml  Iram 
Equation  C-30,  taking  Into  account  that  all  displacements  oro  /oto 
except  ■  V|j.tndujj  *  Vjj,  Note  that  only  those  components 
acting  on  mass  point  22  ora  compntml. 

"  ?rrnlrrrji 1(1 "  ,(v»  ■ v,al  + u  u’811 


5TT+ 


1,000 


Hi 


111  -  0,261(0)  -  0.20(0,0126)) 


•  -2,6klpi 

^*y,ti  “  JfTTT]  “  u'»  *  v»*  "  vi»> 

•  »  •  0  +  0  •  M1MI 

"  -2,Bklpi 

*  anrgRrnra 1,1  ‘  *,><u**  * u,i>  *  "  v,»)1 

•  1,1  ■ M81101  •  0  «»0’aBU 

■  -2,6  klpt 

a.yl,,  *  lu»  •  U|>  *  *l»  - 

■  10  ■  0'0,JI1  ♦ 0  •  01 

*  -2,6  kips 


Nolo  that  lor  tha  first  cycle,  mass  point  22  has  not  heap  moved, 
Hance,  uaJ  *  vja  -  0  ami  all  forco  component*  at  stress  pc'ntw 
21  and  22  ■  0,  Front  syinmolry,  It  can  also  be  cunclti,  ml  Hint 


Stop 


Oporotlon 


The  equality  of  the  shearing  forces  and  the  norma!  forces  at  n  streM 
point  on  this  first  cycle  Is  a  c:o  riddance, 

G  Following  the  sign  convention  of  Figure  C-2 


-(F) 

+  IF.y)„  +  (M„ 

M 

-(-2,6)  - 

(-2,6)  - 

1  +6,0  kips 

“<Mil  “ 

+  fy),,  +  ^sy)|j 

y 

-  (-2,0)  - 

(-2,6)  • 

1  +6,0  kips 

t  Nnwuij  ■  oldUn  +  (U„1W»)|, 

u„  •  0  +  (0,000028X0)  -  3,128  k  10*a  lnoh»» 

Similarly 

v„  -  0  +  (0.000828X8)  -  3,128  m  10*b  InohM 

Not*  that  thaw  displacements  of  mass  point  22  dast1  oy  the 
equilibrium  of  mass  point  12, 

0  This  Is  the  last  mass  point  and  the  and  of  the  first  cycle  of 
relaxation,  Go  to  10, 

10  All  mass  points  are  not  In  equilibrium  since  the  dlsplacam®  its 

U|)  and  V|j  under  step  1  above  destroyed  the  equilibrium  of  mass 
point  12,  Therefore,  there  must  I a  second  cycle  of  relaxation, 
beginning  ,»t  step  d,  Note,  however,  that  In  only  one  cycle  of 
relaxation  and  displacement,  component!!  h„va  outlined  nearly 
00%  of  their  final  velum, 

The  operations  listed  above  demonstrate  the  procedure  for  elastic 
behav!or.  Suppose  that  a  iHifflrlan*  numhor  ol  relaxation  cycles  have  heun 
performed  to  bring  both  must  polmt,  to  within  an  poor  ptablo  error  In  tha 
equilibrium  equations,  The  following  discussion  Indicates  how  the  yield 
criterion  Is  applied  (Step  1 1  of  Figure  C  <11  anil  how  the  force  components 
at  u  yielded  stress  point  (Skip  he  ol  Flyura  C-d)  tire  computed, 


To  illustrate  the  aDDlication  of  the  yle'd  criterion,  assume  that  tho 
yield  stress  in  simple  tension  for  the  material  is  35  ksi.  Then  the  yield  stress 
In  simple  shear  Is 


k  .  .  36  -  17,6  ksl 

k  2  2 

Using  the  equilibrium  displotomonts  from  Equation  C-22  and  applying  the 
yield  criterion  (Equation  C-31)  to  stress  point  12,  wo  get 

-  ^  /kfiV 

(Ujj  -  uia  -  vaJ  +  vn)a  +  (uaa  -  u1#  +  v33  -  via)  <  \-ffj 


(0  -  0.01429  -  0.003571  +  0)3 


+ 


(0,003571  -  0  +  0  -  0.Q1429)3  < 


'(17.BH1.414)  3 

_ 1,000  __ 

2(1  +  0.26) 


0,000434  <  0.003B3 
end  for  stress  point  22  wo  get 

(u33  -  uw  -  v3J  +  vJ3)5  +  (uaa  *  +  vjs  '  v«^  * 

(0  -  0,003671  -  0  +  0}3 


Off 


+  (0  -  0  +  0  -  0,0003571 )3  < 


117, SHI, 414) 


itrimsr 


0,0000255  <  0,00383 


Therefore,  both  stress  points  have  not  yielded  a<  an  ox  ter  no  1  pressure  o( 
14,14  ksl,  r* if st  yielding  will  tfiku  piano  tit  stress  point  1 1  and  12  ot  an 


external  pressure  of 


42  ksl 


B2 


Note  that  this  value  of  external  stress  is  considerably  greater  than  the  yield 
stress  in  simple  tension  or  compression  of  31  ksi  assumed  for  the  material, 
This  Isa  characteristic  of  failure  or  yielding  In  two  dimensional  stress  systems. 

Until  the  external  load  level  has  reached  42  ksl,  all  forces  and 
displacements  Increase  linearly.  When  this  elastic  limit  has  been  reached,  the 
corresponding  forces  and  displacements  are  42/14.14  times  those  of 
Equations  C-33  and  C-34: 


U1J  “  v12  “ 

4.244  x  1Q‘J  Inches 

Ujj  ■  v23  - 

1.061  x  10‘J  inches 

<F*>,a 

-  -23.34  kips 

(Fy)u 

■  -2,12  kips 

(C-20) 


(F,v)ia  -  -8.37  klpi 

These  values  are  recorded  and  are  used  to  determine  the  total  forces  and 
displacements  for  the  first  load  Increment  above  the  42-ksl  load  level, 

Suppose  now  that  the  external  load  level  Is  Increased  to  49.08  ksl.  As 
a  first  approximation  to  the  final  displacements  at  this  new  loBd  level,  the 
displacements  of  Equation  C-30  are  Increased  In  the  same  ratio  as  thn  loads: 

(0,04244)  -  0,04969  Inch 

(C-37) 

(0,01081)  -  0,01240  Inch 

Note  that  two  seta  of  displacements  are  available:  the  last  set  of  equilibrium 
displacements,  Equation  C-30,  end  the  current  set  of  displacements, 

Equation  C-37  (which  In  general  1$  not  compatible  with  the  condition  of 
equlllhr'um),  These  two  sets  of  displacements  are  necessary  In  order  to 
compute  tne  incremental  plastic  force  components  as  discussed  previously, 
These  Incremental  plastic  force  components  can  ue  computed  from 
F.quatljn  C-32  for  stress  point  12  as  follows: 

lF*)ia  *  A(Auaj  -  Auia)  +  B(Avaa  -  Av13) 


uu  ■  v,a  » 

uaa  -  v28  - 


-C(Auaa  -  AUf3  -  Avaa  -  Avu) 
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(Fy)ia  -  B(Au23  -  Au12)  +  A(Av22  -  Av13) 

+  C(Au22  -  Au 1 3  -  Av23  -  Av12) 

IF" *y ) 1 2  ”  -  C(Au23  -  Au,2)  +  C(Av22  -  Av13) 
+  D(Au22  -  Au13  “  Av23  -  Av12) 

G 


1  )  4  G  +  3  K 
where  M  -  j  < - g - 


1 


k2  62 


(120j^3K  +  _0_  nF  ,  . 

}  3  klj2 


<F»'„  ’ 


D  ‘  t[ 


G  -  !Fxv>, 

k2  62  1 


’] 


G 

K 


2(1  +  v) 

E 

3(1  -  2i>) 


Note  that  these  equations  require  Instantaneous  values  of  the  force  components 
If  the  Incremental  force  components  are  to  be  computed.  For  small  Increments 
In  the  external  loading,  the  Instantaneous  forces  are  very  nearly  equal  to  the 
forces  at  the  last  equilibrium  configuration,  Equation  C-36. 

Thus,  the  Incremental  relationships  by  which  the  Incremental  force 
components  In  the  plastic  region  can  be  computed  may  be  determined.  The 
next  step  requires  an  approximation  of  the  forces  at  the  stress  points  by  adding 
the  Incremental  forces  to  the  last  set  of  forces.  Finally,  the  forces  acting  at  a 
yielded  stress  point  are  computed,  and  the  relaxation  technique  proceeds  as 
before, 
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Appendix  D 

DERIVATION  OF  SELECTED  EQUATIONS 
by  W.  D.  Atkins 


In  this  appendix,  the  equations  necessary  to  apply  the  downgrade 
option  in  PERFPLAS  II  are  derived  and  then  the  equations  used  in  subroutine 
CORRECT  are  presented. 

Downgrade  Option 

The  yield  criterion  equation  Is  the  equation  of  a  circle  with  its  center 


Fxy  ■  shear  force  on  the  xy  plane  at  a  stress  point 


6  ■  diagonal  distance  between  mass  points 

k  ■  yield  stress  In  simple  shear  of  material 

Only  the  linear  elastic  cose  will  be  considered,  That  Is,  It  Is  assumed  that  the 
forces  at  any  stress  point  will  vary  linearly  with  variations  In  applied  loads. 
Further  It  will  be  assumed  that  the  forces  at  a  stress  point  which  has  yielded 
at  the  current  loading  are  computed  using  the  elastic  equations,  Referring 
to  Flgu-e  D-1,  the  subscripts  1 , 2,  and  3  denote  conditions  before,  at,  and 
ufter  yield, 

Mathematically 

••  Fy^a  +  “  6k,  yield  criterion 

[(F,i'Fvs)1  +  4F,,jT  " 8k 
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Figure  D-1,  Soil  yield  criterion, 


What  Is  required  to  be  found  Is  the  load  Increment,  AE,  necessary  to 
change  the  state  of  stress  from  condition  3  to  condition  2.  Setting  the  load 
equal  to  L,  then  L3  -  L,  «  AL,  or  change  in  load  and  L2  -  L.j  -  AE,  or 

load  Increment.  If  there  Is  a  change  In  |(F„  -  Fv)2  +  4  Fxv2 1  then  a 
unit  change  In  load  Is  L  ^ 


Subroutine  CORRECT 


At  the  end  of  an  incremen*  of  plastic  straining,  the  stresses  at  an 
already  yielded  point  may  exceed  the  yield  condition,  since  approximate 
linearized  equations  are  used,  For  any  assumed  increment  of  plastic  strain 
along  a  chosen  direction,  the  force-displacement  relations  given  in  Appendix 
C  will  indicate  changes  in  stress  lying  along  the  tangent  to  the  Mohr  circle. 

Thus,  after  several  successive  plastic  strain  increments,  the  stress 
condition  might  well  be  in  excess  of  the  yield  condition.  To  prevent  this,  a 
correction  must  be  made  at  the  end  of  each  increment  to  bring  the  stresses 
back  to  the  yield  surface  along  the  radius  through  the  stress  state  existing  at 
the  end  of  the  increment.  The  correction  to  ho  applied  Is  derived  In  the 
following  pages. 

Referring  to  Figure  D-2,  the  following  quantities  are  known: 

-  u  °*76  °V2S  ’xv2fi 

"2  "  2  '  FV2  "  2  '  F*V2  "  2 

+  ®v,  _  c,x2  +  av2 

2  "  2 
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It  is  nnccsr-'iry  to  define 


and  since  after  adjustment 

J2  ■  k 


i 

l 

1  then  R,  ■  k 


The  derivation  then  proceeds  as  follows:  From  the  low  of  similar  triangles 
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Appendix  E 


DATA  INPUT  GUIDE 
FOR  PROGRAM  PERFPLAS  II 

by 


N.  Shoemaker 


GENERAL  PROGRAM  NOTES 
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Appendix  H 


PHOTOFLASTIC  STUDY  OF  THE  DISTRIBUTION 
OF  MAXIMUM  SHEAR  STRESS  IN  AN  ELASTIC  MEl'lUM 

bv  C.  L.  Liu 


Photoelartlc  methods  of  determining  the  maximum  shear  stress 
distribution  offer  a  convenient  and  Inexpensive  maann  of  supporting  the 
numerical  results  obtalnod  from  the  computer  progiam  u,id  tho  rroature- 

menu  of  the  field  teit  program.  The  circular  polarlscopr  permits  a 
visualization  of  the  shear  stress  patterns  which  develop  ax  loads  are  applied 
to  a  notched  two  dimensional  modal. 

Test  Sptolman  and  Apparatus 

Tha  photoelastic  msdlum  used  In  thaw  tests  Is  urethans  rubber 
having  a  modu'ue  of  alustlclty  of  BOO  psl  and  a  Poisson's  ratio  of  0,40.  This 
Is  equivalent  to  s  value  of  74  psl  for  tha  modulus  of  elasticity  as  measured 
by  performing  trlaxla!  tests  on  soma  nelacted  clay  sediment  cores  obtained 
from  San  Francisco  Bay. 

Loads  wsra  applied  to  the  urethane  rubber  vie  relatively  rigid 
plexiglass  (Luolte)  forms  which  wire  meahinsd  to  precisely  fit  the  notch Jd 
Bras  of  the  urethane  rubber.  The  materials  ware  bonded  by  an  epoxy  cement, 
Difficulties  In  machining  the  urethane  rubber  end  In  bonding  the  two  materials 
resulted  In  the  development  of  residual  stresses  siptiolelly  noticeable  at  the 
Internal  corners,  Aluminum  templates  ware  found  to  bs  well  suited  tor  making 
the  models,  See  Figure  Ho  for  the  general  specimen  dimensions. 

Tha  test  apparatus  used  In  the  experiment  Is  shown  sohsmitloelly  In 
Figure  H«2.  Briefly,  the  Isochromatlc  pottern  of  •  specimen  strained  In  the 
steal  ingle  frame  Is  viewed  through  the  polerlsQope  analyzer,  and  the  Image 
Is  racordod  bv  the  Slner  camera,  A  simple  mechanical  jack  Is  used  to  position 
the  specimen,  Figure  H<3  Illustrates  how  the  tension  force  on  the  spaclmen  Is 
controlled  by  a  w„*,g  nut,  Tha  tensile  force  on  the  modal  was  animated  by 
reproducing  tha  stress  pattern  with  a  known  foice  applied  bv  a  simple  lever 
balance  system,  as  shown  In  Figure  H-4, 

Procedure 

The  calibration  constant  C  Ik  the  ratio  of  the  moxlmum  shear  stress, 
rnim,  to  tha  order  of  birefringence  R  of  tha  Isochromatlc  pattern,  A  2-Inch 
circular  urathana  rubbar  disc  was  uaed  tor  ndlbrutlon,  A  photograph  of  tha 


disc  loaded  in  pure  compression  by  4,5  pounds  is  shown  in  Figure  H-5.  The 
center  of  the  disc  Is  the  reference  point.  The  birefringence  order  at  this  point 
is  6.0.  The  value  of  rmm  at  the  center  was  calculated  analytically:  rmm  - 
4  F/ir  t  D .  The  compression  force  is  4,5  pounds,  the  model  thickness,  t,  is 
0,25  Inch,  and  the  model  diameter,  D,  is  2  inches.  The  value  of  rmm  at  the 
reference  point  is  found  to  be  1 1 .5  psi.  Thus,  the  value  of  C  is  equal  to 
1.91  psi.  To  aid  In  data  reduction,  th6  linear  shear-birefringence  relation  is 
shown  In  Figure  H-6. 

Photographs  of  unloaded  specimens  revealed  Internal  stresses  which 
were  found  to  exist  near  the  notched  internal  corners.  These  specimens  were 
then  tested  with  loads  sufficient  to  produce  clear,  distinctive  birefringences 
of  a  reasonable  quantity.  A  half-embedded  circular  cylindrical  specimen  was 
subjected  to  approximately  equal  load  Increments  until  failure,  and  sequential 
photographs  were  taken.  The  force  levels  in  each  case  were  estimated  by 
reloading  the  specimen  with  the  lever  balance  system  so  that  the  Isochromatlc 
pattern  In  the  photographs  was  reproduced, 
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Figure  H-1,  Typical  dlmanalona  of  aoll  modak 


f'lgur*  H-3.  8p*olm»n  toadr. 


1,000  g 


Flgurt  H-4.  Foret  mttiurtmtnt  (litoll, 


All  nf  the  photographs  were  taken  with  black  and  white  Polaroid 
Land  type  57  films  and  a  Wrattan  77a  yellow  filter.  It  was  found  that  the 
best  dark-field  polariscope  picture  was  made  by  using  3000  ASA  film  with  an 
f32  opening  and  an  exposure  of  7  seconds.  A  blue  monochromatic  light 
source  was  used  since  it  gave  more  distinctive  black  and  white  photographs. 

In  all  of  the  photographs,  a  transparent  scale  was  used  as  a  length  reference. 

Results 


A  summary  of  model  dimensions  is  shown  in  Figure  H-7,  Figures  H-8 
to  H- 1 5  shows  a  group  of  models  under  tensile  loadings.  Sequential  photo¬ 
graphs  of  model  F  under  stepwise  load  increments  are  shown  in  Figures  H-16 
to  H-18.  Each  loading  Is  given  in  Table  H-1. 

Interpretation  and  Application 

Based  on  the  photoelastic  principle,  the  maximum  shear  stress,  rwm, 
is  proportional  to  the  birefringence  order  N.  The  rmm  distribution  is  clearly 
shown  In  each  of  the  isochromatic  photographs.  Each  black  curve  may  be 
assigned  a  birefringence  order  along  which  the  value  of  rmm  Is  constant 
Qualitative  conclusions  can  be  inferred  directly  from  the  photographs; 
however,  quantitative  conclusions  require  a  determination  of  the  proportional 
or  calibration  constant, 

By  observing  the  growth  of  the  birefringence  pattern  during  the  test, 
one  may  determine  the  direction  of  the  birefringence  variation,  Referring  to 
Figure  H-9,  the  curved  arrows  Indicate  the  directions  of  increasing  birefringence. 
For  example,  the  value  of  Tmm  Is  zero  at  the  free  corner,  is  low  at  the  center  of 
the  Interface,  and  is  highest  at  the  lower  corners  of  the  Joint,  Consequently, 
failure  Is  Initiated  at  the  Internal  corners  and  would  extend  downward  and 
inward  following  the  "ridge"  of  the  rmm  contour  to  form  a  failure  arc,  which 
Is  Indicated  with  a  dashed  line, 

Failure  curves  for  the  other  boundary  geometries  ere  also  indicated  in 
the  pertinent  photographs.  As  shown  in  the  sequential  photographs, 

Figures  H-16  to  H-18,  the  failure  curve  is  Independent  of  t  ie  external  loading 
for  a  given  boundary  geometry,  Behavior  in  the  plastic  range  cannot  be 
simulated  with  the  urethane  rubber  models, 

In  order  to  make  some  quantitative  estimates  of  maximum  shear  stress 
under  a  given  loading,  the  following  formule  Is  useful' 

rmm  -  CN  (H-1) 
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where  rmm  is  the  maximum  shear  stress  in  the  model  at  a  point  coi responding 
to  tho  birefringence  order  N  and  C  is  the  calibration  constant.  Thus,  for  any 
point  in  the  medium,  the  maximum  shear  st  ess  may  be  calculates.  To 
extrapolate  the  results  to  field  conditions,  the  following  formula  Is  given: 


T  rr 

rm  ■- 1 


(H-2) 


where  m  =  model 

F  ■  applied  force 
L  -  length 
t  ■  thickness 

r  ■  maximum  shear  stress  In  prototype 

In  order  to  compute  the  breakout  force,  rearrange  Equation  H-2  into 
the  following  form: 


F 


(H-3) 


The  procedure  requires  ( 1 )  selection  of  a  critical  point  In  the  medium  and 
(2)  determination  of  the  level  of  loading  Fm  at  which  the  maximum  shear 
stress  reaches  the  yield  stress  at  the  critical  point.  The  critical  point  Is 
located  at  the  center  of  the  failure  ridge.  The  breakout  force  Is  easily 
calculated  from  Equation  H-3, 

As  an  example,  suppose  It  is  desired  to  find  the  mexlmum  shear  stress 
at  a  point  36  Inches  left  of  the  centerline  and  24  Inches  below  the  sediment 
rurfacc  of  a4x4x  16-foot  parallelepiped  which  Is  embedded  6  Inches  In  the 
soil  under  an  applied  load  of  10,000  pounds.  This  embedment  condition 
corresponds  to  the  boundary'  geometry  of  Figure  H-0,  where  the  equivalent 
point  is  designated  as  point  A,  The  birefringence  order  N  at  A  Is  4  6.  From 
Figure  H-6,  rmm  -  8,6  psl,  and  from  Table  H-1 ,  F  •  6,14  pounds,  so  that 
from  Equation  H-2 


10, 000(0, 751(0.261(8.0) 
r  "  (6.1 4)(48) ( 1  2) 


0.303  pel 


which  is  the  predicted  maximum  shear  ntrees.  As  another  example,  suppose 
we  are  interested  In  determining  the  breakout  force  of  a  4x  4x  16-foot 
parallelepiped  embedded  6  Inches  in  a  soli  having  a  shear  yielding  stress  of 
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0,6  psl.  First,  th8  critical  point,  which  is  point  B  in  Figure  H-9,  must  be 
found,  With  a  birefringence  order  of  6,  the  value  of  r„„  is  found  from 
Figure  H-6  to  be  equal  to  11.5  psi,  Using  Equation  H-3,  we  have 


0,6(6, 14)(48)(150)(12) 
<0.75)(0.26)(1 1,6) 


14,750  pounds 


This  value  will  be  higher  than  the  actual  force  since  failure  does  not 
occur  simultaneously  at  all  points  along  the  failure  arc,  During  the  Failure 
process,  the  transfer  of  load  along  the  failure  ridge  causes  the  stress  at  the 
critical  point  to  attain  the  yield  point  much  earlier  than  predicted,  However, 
the  omission  of  the  plastic  consideration  does  not  detract  from  the  value  of 
the  study  reported  on  herein.  It  is  useful  to  ( 1 )  support  the  numerical 
technique  given  in  Appendixes  C  through  G,  (2)  visualize  the  mode  of  failure, 
and  (3)  present  approximate  engineering  estimates  of  the  breakout  force. 


Table  H-1 .  Calculation  of  Forces  on  Models 
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Floui't  H-7.  Gaomitrlii of  modali. 
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Figure  HOI.  Modal  D,  vertically  loaded. 


Figure  H02.  Model  E,  vertically  loaded 


lure  H-0.  Modal  B,  vertically  loaded. 


Figure  H*  10.  Model  C,  vertically  loaded, 
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eiQuri  H<16.  Model  P,  S  2Q'pound  tonilon. 
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